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Im Gegensatz zu den Versuchen anderer Autoren, das 
Neumannsche Randwertproblem für das Rotations- 
ellipsoid mittels der Lam&-Funktionen zu lösen, erhält 
man hier eine Lösung dadurch, daß das Rotationsellip- 
soid auf die Soldnersche Bildkugel abgebildet wird. Bei 
den geodätischen Anwendungen genügt es, in den Ab- 
bildungsgleichungen nur die in der Erdabplattung line- 
aren Glieder zu berücksichtigen. Die Lösung erzielt man 
in Form von Korrekturen, die an dem Ergebnis der sphä- 
rischen Problemstellung anzubringen sind. 


In der zweiten Arbeit wird eine mathematische Lösung 
für das Problem gefunden, aus einer Beobachtung eines 
nahen Satelliten die Lage einer „Standlinie“* für die 
Position des irdischen Beobachters anzugeben. 


In der dritten Arbeit wird gezeigt, wie sich die Elemente 
des Maßtensors eines in der Geodäsie vorgeschlagenen 
nichteuklidischen Koordinatensystems durch geodäti- 
sche Messungen ermitteln lassen. 


Die vierte Arbeit gibt für die im Titel genannten Größen 
geeignete mathematische Entwicklungen an. 
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Spannungsoptische Untersuchungen 
schiefwinkliger elastischer Platten”) 


Von H. SCHWIEGER und G. HABERLAND**) 


verfahrens die maximalen Torsionsmomente und die Richtungen der Hauptbiegungsmomente bei einer 
schiefwinkligen Platte in Form eines gleichseitigen Parallelogramms bestimmt. Die Parallelogrammplatte 
‚konstanter Dicke ist nur an zwei gegemüberliegenden Rändern frei drehbar gestützt und in der Mitte durch 
eine Einzelkraft belastet. Aus den spannungsoptisch erhaltenen maximalen Torsionsmomenten und den 
obigen Richtungswinkeln werden schließlich nach einem geeigneten Integrationsverfahren die einzelnen 
Biegungsmomente errechnet. E 


It is shown how, for a plate having the form of an equilateral parallelogram, the maximum twisting mo- 
ments and the directions of the prineipial bending moments may be determined by employing the photoelastic 
two-sheet method. The plate, which is supposed to be plane-parallel, is simply supported at two opposite 
edges, and has a single load applied at the centre. From the above-mentioned moments and directions deter- 
mined by the photoelastic method, the principal bending moments are finally determined, using a suitable 
integration method. 


B 3To0ü cTaTpe IIOKA3aHO, KAaR C IIOMOINBIO OITUHYECKOTO METONaA ABYX CJIOEB MOSKHO olpe- 


f HeJIHTbB MARKCHMAJIBHOE MOMEHTEBI KPy4YEeHUA H HAaupPpaBJIeEHHUA TJIABHbIX U3THÖAIINHX MOMEHTOB 


IIIACTUHRH, HMeromeii PopMy PoM6a. IlsIactuHKa HOCTOAHHOU TOJIIHHBI CBOOGOAHO OUHPAeTCH 
Ha ABA IPOTHBOIIOJIO}KHEIX KPan, B CepelIuHHe HAXONHTCAH HAarpy3Ka. Vlcxo]A H3 HOJIYYeHHBIX 
OITHYECKUM METOIOM MAKCHMAJIbHbIX MOMEHTOB KPy4YeHHA HM BbIlIe YIIOMAHYTBIX YTJIOB HA- 
IPaBJIeHUA, UHTETPUPOBAHHEeM BbIYHCHIAIOTCH H3TMHÖAIOIIME MOMEHTLI. 


1. Einleitung 


Schiefwinklige, elastische Platten sind in neuerer Zeit häufig im Brückenbau, wo sie als 
Fahrbahnplatten dienen, verwendet worden. Entsprechend dieser technischen Bedeutung sind 
die Probleme an schiefwinkligen Platten in zunehmendem Maße Gegenstand sowohl theoretischer 
wie auch experimenteller Untersuchungen. Das umfangreiche Schrifttum ([1] bis [36]) bezeugt 
deutlich, welches Interesse man den parallelogrammförmigen oder schiefwinkligen Biegeplatten 
entgegenbringt. Bei der theoretischen Behandlung bedient man sich vorwiegend der Differenzen- 
rechnung zur angenäherten Berechnung der verschiedenen Plattengrößen, wie Durchbiegung 
und Biegungsmomente. Auch hat man die partielle Differentialgleichung für die Plattendurch- 
biegung und die einzelnen Plattengrößen in schiefwinkligen Koordinaten dargestellt. In diesem 
Falle hat die Plattendifferentialgleichung durch zusätzliche Glieder einen komplizierteren Auf- 
bau gegenüber der Darstellung in kartesischen Koordinaten, jedoch gewinnt man eine bessere 
Übersichtlichkeit bei den Randbedingungen. Bei den experimentellen Methoden zur Unter- 
suchung elastischer Platten überhaupt bedient man sich neben den spannungsoptischen Ver- 
fahren noch der Krümmungsmessungen und Dehnungsmessungen. Auch finden Spiegelanordnun- 
gen und das Moire-Verfahren Anwendung zur Analyse des Spannungszustandes der Platte 


- [(S. 37] bis [53]). 


Im folgenden wird nun eine Schilderung der Versuche an schiefwinkligen Platten von der 
Form eines gleichseitigen Parallelogramms (Plattenwinkel & = 60°) gegeben. Die Parallelo- 
grammplatte war dabei an zwei gegenüberliegenden Rändern frei drehbar gestützt, an den 
anderen gegenüberliegenden Rändern ungestützt und in der Mitte durch eine Einzelkraft be- 
lastet. Zur experimentellen Untersuchung benutzten die Verf. das spannungsoptische Zwei- 
schichtenverfahren, welches schon früher von den Verf. bei anderen Plattenproblemen erfolgreich 
angewendet wurde. Auch wird aufgezeigt, wie durch ein der Plattenform angepaßtes Auswerte- 
verfahren unter Benutzung der spannungsoptisch bestimmbaren maximalen Torsionsmomente 

*) Vorgetragen auf der GAMM-Tagung in Freiberg/Sa., 19.—23. April 1960. 

**) Institut für Angewandte Mathematik und Mechanik der DAW zu Berlin, Spannungsoptisches 


Laboratorium. 
19 


In der folgenden Abhandlung wird geschildert, wie man mittels des spannungsoptischen Zweischichten- 


z 


a u Ele 


a \ 


und der Richtungen der F Jiegur 
erhalten wurden. Die Ergebnisse der sehr v 
der verschiedenen Biegungsmomente, wurden in Fo 
Platte festgehalten. 


i \ Br 32 5 Baar: , 
re Pe) BEREIT 
2. Elastooptisches Verhalten von Biegeplatten — 

Die spannungsoptische Methode geht in ihrem Wesen auf die Spannı 
zurück, wonach viele durchsichtige, von vornherein optisch isotrope Stoffe, wie 
Kunststoffe, bei mechanischer Beanspruchung optisch doppelbrechend werde 
nungsoptischen Untersuchung von elastischen Platten muß man allerdings | | 
kehrungen treffen, um einen optischen Effekt zu erhalten. Durchstrahlt man näm 
Biegeplatte, die innerhalb der KırchHorrschen Bedingungen belastet ist, in Richtung ihr. 
malen mit linear polarisiertem Licht, so bleibt bei Beobachtung dieser Platte durch de 
Polarisator gekreuzten Analysator das Gesichtsfeld dunkel. Das einfallende linear polarisi te 
Licht wird zwar beim Eintritt in die Platte wegen der Spannungsdoppelbrechung in zwei ib > 
gonale Lichtkomponenten aufgespalten, und diese haben auch eine optische Phasendifferenz 
beim Durchgang bis zur Plattenmittelfläche miteinander, jedoch wird die optische Phasendiffe- 
renz beim Durchgang durch die andere Plattenhälfte wieder bis auf Null ‚vermindert. Wegen 
des Vorzeichenwechsels der Biegungsspannungen an der Mittelfläche sind die optischen Phasen- 
differenzen der beiden Plattenhälften entgegengesetzt gleich, so daß sie sich für die ganze Platten- 
stärke gegenseitig kompensieren. Bei weiterer Durchbiegung der Platte kann man allerdings 
einen spannungsoptischen Effekt beobachten, da sich dann dem Biegungsspannungszustand ein 
Membranspannungszustand überlagert. Die spannungsoptische Beobachtung gibt also eine gute 
Möglichkeit zu entscheiden, inwieweit, gekennzeichnet durch gewisse Aufhellungen des Gesichts- 
feldes, die Platte reines Biegungsverhalten zeigt, also die KırcaHorrschen Plattenbedingungen 
noch erfüllt sind. Die Aufspaltung des linear polarisierten Lichtes erfolgt, wie oben schon er- 
wähnt, in zwei linear und senkrecht polarisierte Lichtkomponenten, deren Schwingungsrichtun- 
gen mit den Richtungen der Hauptbiegungsmomente zusammenfallen. Diese Richtungen 
werden auch über die Plattenhöhe beibehalten, da die sekundären Hauptspannungen?!) nur 
durch die dem Momententensor zugehörigen Spannungen gegeben sind. Die sekundären Haupt- 
spannungen fallen also mit den Hauptspannungen des Biegungszustandes zusammen. (Die 
Schubspannungen der Querkräfte haben keinen Einfluß auf Schwingungsrichtung und optische 
Phasendifferenz der beiden Lichtkomponenten.) 


Zur Erzielung eines spannungsoptischen Effektes setzt man nun die Modellplatte aus zwei 
Schichten durchsichtiger Werkstoffe mit unterschiedlicher optischer Wirksamkeit zusammen 
(Zweischichtenverfahren [54] bis [59]). Auf diese Weise können sich die optischen Teileffekte 
in den beiden Schichten nicht mehr aufheben, und man kann durch den Analysator Isochromaten 
und Isoklinen beobachten. Dabei sind die Isochromaten, welche bei Verwendung einfarbigen 
Lichtes als dunkle Linien erscheinen, Orte gleicher optischer Phasendifferenz. Im weißen, 
linear polarisierten Licht sind die Isochromaten Orte gleicher Interferenzfarbe, so daß ent- 
sprechend der spektralen Zusammensetzung des weißen Lichtes die einzelnen Streifen der ver- 
schiedenen Ordnungen farbig aufgegliedert sind. Beim Drehen des Polarisationskreuzes ändern 
die Isochromaten ihre Lage nicht. Die Isoklinen, welche sich sowohl im einfarbigen wie mehr- 
farbigen Licht stets als schwarze Linien abzeichnen, verschieben sich dagegen bei inhomogenen 
Spannungszuständen im Gesichtsfeld beim Drehen des Polarisationskreuzes. Längs der Iso- 
klinen fallen das eingestellte Polarisationskreuz und das Kreuz der Hauptspannungsrichtungen 
zusammen, so daß also keine Komponentenzerlegung des einfallenden linear polarisierten 
Lichtes stattfindet und demzufolge der Polarisationszustand des polarisierten Lichtes nach 
dem Durchgang durch die belastete Platte unverändert bleibt. Der Isoklinenparameter — der 
Richtungswinkel eines Hauptbiegungsmomentes @ (bzw. @ + 90°) — wird dabei auf das ge- 
wählte x, y-Koordinatensystem bezogen, in welchem dann 9 = 0° bzw. 90° ist. 

Wie sich ableiten läßt, sind bei Platten mit konstanter Dicke die Isochromaten die Orte 
gleichen maximalen Torsionsmomentes. Zwischen der auf ?x bezogenen optischen Phasen- 
differenz (relativer optischer Gangunterschied) m und dem maximalen Torsionsmoment T 
besteht folgender Zusammenhang: 


max 


m Tax , 48 n (K,—K}) 


h (m-+l4n+ 1) 
Diese Beziehung gilt für den Fall zweier gleich starker Schichten (h/2) und zeigt deutlich, 
daß bei gleichem elastooptischem Verhalten der Schichtmaterialien — gekennzeichnet durch 


= 0. (MM) 0 en 


!) Siehe M. M. Frocur, Photoelasticity, Vol. II, p. 333, J. Wiley and Sons, New York 1948. 
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die Gleichheit der spannungsoptischen Konstanten K, und K, — der optische Effekt verschwindet. 
Der Einfluß der Elastizitätsmoduln der Schichten wird durch das Verhältnis n = E,/E, aus- 
gedrückt. Weiterhin ist Gleichung (1) für den Fall gleicher Poıssoxscher Konstanten der be- 
nutzten Modellwerkstoffe abgeleitet worden. Weichen die beiden Poıssoxnschen Konstanten 
voneinander ab, so sind gesonderte Betrachtungen?) notwendig. 


Da man bei Platten nach dem spannungsoptischen Zweischichtenverfahren im allgemeinen 
nur wenige, relativ weit auseinander liegende Isochromaten bekommt, die allerdings durch 
Veränderung der Belastung in gewissen Bereichen kontinuierlich verschoben werden können, 
so ist es doch im allgemeinen erforderlich, kleine m-Werte in den weniger beanspruchten Platten- 
gebieten nach den bekannten Kompensationsmethoden zu bestimmen. Bei den nachstehend 
wiedergegebenen Versuchen wurden die kleinen optischen Gangunterschiede nach der Methode 
von TArpy [60] mittels A/4-Platten gemessen. Zum besseren Verständnis des eben Aufgeführten 
und auch über die Grundlagen der Spannungsoptik sei auf die allgemeine Literatur am Schluß 
dieser Arbeit verwiesen. 


3. Das Auswerteverfahren 


Wie vorne ausgeführt, kann man aus den spannungsoptischen Bildern sowohl das maximale 
Torsionsmoment Tax wie auch die Richtungswinkel bzw. @ + 90° der Hauptbiegungs- 
momente M, bzw. M, entnehmen. Aus diesen spannungsoptischen Ausgangsgrößen läßt sich 
sofort das Torsionsmoment 


Vay > mas SID Age Pe er (>) 


berechnen. Diese Plattengröße, insbesondere deren Abhängigkeit vom Ort, bildet den Ausgang 
eines Auswerteverfahrens [61] bis [63] zur separaten Bestimmung der einzelnen Hauptbiegungs- 
momente M, und M, bzw. der Biegungsmomente M, und M,. Erinnert man sich, daß die 
einzelnen auf 1 cm Schnittlänge bezogenen Momente mit den örtlichen Ableitungen der Durch- 
biegung der Plattenmittelfläche w(z, y) in folgendem Zusammenhang stehen: 


o2w o:w 

ri | at en Ba 92, 
ww O2 

ag (D tu | AN 
2w 

Ty=-B( nn ar Sa aa (5) 


(x, y-Ebene fällt mit Plattenmittelfläche zusammen), dann erhält man zunächst durch Inte- 
gration der Gleichung (5) in x- bzw. y-Richtung 


> 
2, 


ow 1 'ow 
a el er ee IT BORN SSR REN 5 0 2 PRRREREN EN 
Fun B (il =» j "4 ; (), 
Po 
P, 
dw 1 ow‘ 
ee | Ted = a ET 
&).. Bi—u) ih wäyT (ür)r. @) 
p 
und nach anschließender Differentiation 
P, 
Fw N oT dx - (3) RT PERS), 
(ar P, BI —-u) 0 Oy? >, 
P) 
P’ 
a, 2 (=) er SAN 0). 
(ae), Bi —u) Gi 0x Ps 


2) B. Gıre [56]. 
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w ow 
Führt man hier die Ersatzmomente m, =— B(1—.ı%) Fr undm, = — B (18) Er 
.. ein, so ergibt sich: 
P, 
(m)r,=(+u) on m FÜ), a Ba 
Po . 
P, 
1 Par, d a 
(m,)r, = ( 4 | ara + tm SRarE Ir I 2 n 


Pr 
In den vorstehenden Gleichungen bedeutet B=E RP/12(1—#) dieBiegungssteifigkeit, welche 
sich aus dem Elastizitätsmodul E, der Poıssonschen Konstanten « und der Plattenhöhe h er- 
rechnet. Bei einer Zweischichtenplatte ergibt sich eine Biegungssteifigkeit, die sich aus den 
elastischen Konstanten beider Schichten und deren Dicken zusammensetzt (S. [56]). Wie man 
den Gleichungen (10) und (11) entnimmt, benötigt man 
für die Auswertung nur die Poıssonsche Konstante. 


03cm 


Im Falle der schiefwinkligen Platte ist es nun 
zweckmäßiger, nicht die Differentialquotienten OT’, „/oy 
zu bilden, denn diese lassen sich an den Plattenrändern 
schlecht ermitteln, sondern die Differentiation in u-Rich- 
tung vorzunehmen. Die Lage des Koordinatensystems 
und auch ein willkürlich gewählter, geradliniger Inte- 
grationsweg sind in Bild 1 wiedergegeben. Mit den Be- 
zeichnungen dieses Bildes bildet man 


Tas. OD Mi: RE ANETTE 
du 0x du 0%y du WERE °y na 


Bild 1. Schema der zweiseitig freidrehbar aufgelagerten (12), 
Modellplatte mit einer Schiefe von « = 60°. P,ist ein und 
Ausgangspunkt und P, ein Endpunkt bei dem Inte- 


grationsverfahren. Das Koordinatensystem x, y liegt oT 1 dT,, EP> 
in der Plattenmittelfläche = — cotga AR 
oy sina du 0x 


Durch Einsetzen in Gleichung (10) bekommt man schließlich 


2, 


| Ay ® 2; in | ein de — {(Tz Pr, a] (Typs Ric ( 


m, 
nt A 
1+ -E (14) 


ru 


l+u - sin & du 
Po 


An den Ausgangspunkten der Integration P,, die man an den Plattenrand legt, lassen 
sich T,, und m, aus den spannungsoptisch bestimmten Werten von Tax berechnen und dem- 
zufolge Teile der Gleichung (14) zusammenfassen’). Für Vorzeichenüberlegungen ist es aller- 
dings besser, die Form der Gleichung (14) zu benutzen. 

1 


. m h 
Aus dem Ersatzmoment es lassen sich dann einmal 


ne 37 = tr 2 Tax 608290: 27,202 ai (15) 
und andererseits weiterhin 
> 

M,= 7 (my BB) 5 IR En (16), 

IM i 

M,= II® (my +. ma): 1a. Se 
berechnen, 

®) Es gilt am freien Rande ) = “ (Imaz)pı 
P, 


HM 


37 (sin?« — u. cos®«) und (77y)p, = (Tax) * sin 2a. 
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4. Versuehsdurchführung und Ergebnisse 


Zur Herstellung der Zweischichtenplatten wurden einmal dünne Plexiglasplatten und zum 
anderen Platten aus dem tschechoslowakischen Modellwerkstoff Umapolar (BZM) bzw. aus 
dem Werkstoff VP 1527 (Dynamit AG., Köln-Troisdorf) benutzt. Diese Platten wurden sowohl 
nach dem Gesichtspunkt großer Freiheit 
von Eigenspannungen wie auch nach ihrer 
höchstmöglichen Planparallelität ausge- 
sucht und dann unter leichtem Druck mit 
dem Polymerisationskleber Piacryl SH 
(Stickstoffwerk Piesteritz) zusammenge- 
klebt. Vermutlich durch die bei der 
Polymerisation frei werdende Wärme ent- 
stehen in den Zweischichtenplatten polari- 
sationsoptisch beobachtbare unregelmä- 
Bige Aufhellungen (,‚Wölkchenstruktur‘‘), 
die aber im Laufe der Zeit mehr und 
mehr verschwinden. Die nahezu gleich- 
seitigen, schiefwinkligen Modellplatten 
mit einem Plattenwinkel von 60° hatten 
nach der mechanischen Bearbeitung fol- 
gende Abmessungen (vgl. dazu Bild 1): 


Modellplatte aus Plexiglas-BZM: 


a=12 (em), ‚h : a Bild 2. Belastungsanordnung 
Modellp latte aus Plexiglas-V = A: Analysator mit A/4-Platte, Pl: Modellplatte, St: Stützlager, 
a = 16 (cm), h= 1,09 (cm). F: Ringdynamometer mit Druckstempel 


Pa SAhn) 


a) P=50(kp), b) P = 100 (kp) 
4 2 BIeslark] atte ‚stehe aus xielas aus de schechoslowakischen Werkstoff Uma- 
3. Isochromatenaufnahmen der Zweischichtenplatte, bestehend aus Plexiglas und aus dem tschechoslowa 1 
lar (BZM) bei zwei Laststufen zwischen gekreuzten Polarisationsfiltern (Zirkularfeld). Die Zahlen geben die Isochromatenordnungen an. 
Die Einzellast wirkte in Plattenmitte auf die BZM-Schicht. Der untere Rand ist durch die Auflagerrolle verdeckt (A = 589 (ma)) 


b)e = 15° 
; i - 3 > »s Polarisationskreuzes. (9 = 0° bzw. 90° fällt mit der x- bzw. y-Achse zusammen.) 
. , Isoklinenaufnahmen bei verschiedenen Stellungen des Polarisationskreuzes. (g % bzw. 1 it der w. y-A 
nn Die Zweischichtenplatte besteht in diesem Fall aus Plexiglas und VP_1527 (Dynamit-AG. Köln-Troisdorf). 
Die Einzellast wirkte auf die Plexiglasseite 


a)9= 0), 
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Zwecks Auflagerung wurden zwei Parallelseiten etwas länger bemessen, damit zv on 
den Stützrändern und den freien Rändern die Platte die Form eines gleichseitigen Fi lo- 
gramms hatte. Ein Stützrand war als Schneide ausgebildet und der andere war ein Rollenlager. 
Die Einzelheiten der Belastungsanordnung und auch ein Teil der spannungsoptischen Bank sind 
Bild 2 zu entnehmen. Die Belastungsanordnung befindet sich dabei zwischen drehbar gekoppel- 
ten, gekreuzten Großflächen-Polarisatoren, und die Ausleuchtung des Gesichtsfeldes erfolgt 
sowohl mit weißem wie auch’ mit monochromatischem, diffusem Licht. In den Bildern 3 und 4 


sind als Beispiele einige Isochromaten- und Isoklinenaufnahmen festgehalten. Die Aufnahmen 


(7 217 
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Bild 5. Zusammenfassende Darstellung der Isoklinen. Bild 6. Trajektorien der Hauptbiegungsmomente. 

Der Verlauf der Isoklinen an den Stützrändern ist teil- Die verschiedenen singulären Punkte sind mit $, bis 5, 

weise gestrichelt gezeichnet, da der genaue Verlauf nicht bezeichnet 
ermittelt werden konnte 


wurden mit einer langbrennweitigen Plattenkamera aus größerer Entfernung von der Modell- | 
platte (wegen der Kriecherscheinungen nach einer Belastungsdauer von ca. 10 Min.) aufgenom- 

men. Eine zusammenfassende Darstellung der Isoklinen ist in Bild 5 und die daraus gezeichneten 
Hauptmomentenlinien sind in Bild 6 wiedergegeben. Am Rande ist ein Teil der Linien ge- 
strichelt gezeichnet, da sich in diesen Gebieten wegen des optischen Randeffektes und der opti- ° 
schen Wirkung der Auflagerkräfte die Iso- 
klinen nicht genau ermitteln lassen. 

Zur Festlegung der Isochromaten 
als Tmax-Linien ist ein Eichversuch nötig, 
der daraufhinausläuft, die wellenlängenab- 
Bild 7. Isochromaten des zweischichtigen Eichstabes (Plexiglas — Bzw). Nängige Materialgröße € der Gleichung (1) 
Der Balken N a e a durch eine zu bestimmen. Dazu wurden aus den zu- 
sammengeklebten Platten nicht nur die 
Modellplatten, sondern auch Eichbalken 
mit einer Breite von d = 1 (cm) heraus- 
geschnitten. Diese wurden dann auf zwei 
Stützen gelagert und in der Mitte durch 
eine Einzellast P belastet. In diesem Falle 
stieg das Biegungsmoment M, vom Auf- 
lager zur Laststelle hin linear an (M, ist 
gleich Null), so daß sich in äquidistenten 
Abständen Ax zueinander parallele, senk- 
recht zur Balkenachse verlaufende Iso- 
chromaten ergeben, wie sie in Bild 7 fest- 
gehalten sind. Bei dem Eichversuch be- 
Bild 8. Linien gleichen Torsionsmomentes Tax in (kp) für eine Einzellast fand sich die Biegeeinrichtung ebenfalls 

von P=1 (kp) in Plattenmitte zwischen gekreuzten Polarisationsfiltern 
und wurde parallel zur Einzellast mit 
linear polarisiertem Licht durchstrahlt. Von einer Isochromatenordnung zur nächsten ist 
Am =1. Die Konstante C wurde schließlich aus der Beziehung 


cd 2d Am 
Jo P % Ax = = . ” . . . . * . . * . . ” (18) 
berechnet und erbrachte für die einzelnen Schichtenkombinationen bei Na-Licht (1 = 589 (mu)): 
Balken Plexiglas-BZM: C = 0,134 (Ordn./kp) 


Balken Plexiglas-VP 1527: C = 0,0617 Ordn./kp). 


Nunmehr war es möglich, die verschiedenen Isochromatenaufnahmen auszuwerten. Zur Aus- 
SErUnS wurden beide Zweischichtenplatten herangezogen und von den spannungsoptischen 
Werten Mittelwerte gebildet und alle durch den Zahlenwert der jeweiligen Einzellast dividiert, 
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so daß die Ergebnisse für die Belastung von 1 (kp) gültig sind. Das Ergebnis für das maximale 
Torsionsmoment Tax ist in Bild 8 festgehalten. Anschließend wurden die Torsionsmomente 
T,, berechnet und für hinreichend viele Schnitte in u-Richtung aufgetragen und differenziert. 
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= = const für P=1(kp) Bild 10. Kurven 7 
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ar = const für P=1(kp) 
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Bild 13. Linien gleicher Hauptbiegungsmomente M, 


Durch Integration nach Gl. (14) wurden für verschiedene Wege parallel zur x-Achse die Ersatz- 
momente m, bestimmt unter Berücksichtigung der Randwerte in den jeweiligen Punkten P, 
(S. Bild 9). Bei der Auswertung war es notwendig, den Vorzeichen der einzelnen Summanden 
der Gleichung (14) besonderes Augenmerk zu schenken. Die Ergebnisse der umfangreichen Aus- 
wertung findet man in den Bildern 10—16, wobei aus Timax, 9 und m, alle anderen Platten- 
momente nach den Zusammenhängen des Mourschen Kreises berechnet wurden. 
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5. Schlußbemerkungen 


Eine Kontrolle über die erhaltenen Werte der Plattenmomente ist möglich, wenn man 
einmal die Integration nach Gleichung (14) vom rechten bis zum linken Plattenrand durchführt 
und zum anderen das Gleichgewicht an einem endlichen Plattenstreifen ausrechnet. Im ersten 
Falle muß sich wieder an der oberen Grenze der Integration ein Ersatzmoment ergeben, welches 
sofort aus dem spannungsoptisch bestimmbaren Tyax-Werten am Rande berechnet werden 
kann. In diesem Sinne kann man Gleichung (14) umstellen und erhält mit einer Poıssonschen 
Konstanten von u = 0,35: 

Pa 
1,28 - {(Timax)p_ — (Tmax)Pp } = Toy gr BE II (LI 
’ ax, a aAX f) du 


Pr, 


(Im folgenden liegt P, am linken und P, am rechten Plattenrand.) 
Die Ausführung der Rechnung für eine Platte mit a = 12 (cm) erbrachte die in Tabelle 1 


wiedergegebene Gegenüberstellung. Wie man dieser Tabelle entnehmen kann, ist die Überein- 
stimmung recht gut (Spalte 4 und 5). 


Tabelle 1 
ß er 2 Pa 
| | 1,28 S(Tmax)P f dTz 
\p L Ey dx 
77 (T max) D, (T max)? 2 ei Tas.) ie da 
| P 
6,9 (cm) | 0,056 (kp) \ 0,056 (kp) 0 (kp) 0 (kp) 
8 (cm) | 0,045 (kp) | 0,067 (kp) 0,028 (kp) 0,028 (kp) 
9 (em) | 0,033 (kp) | 0,083 (kp) 0,064 (kp) 0,065 (kp) 
10 (cm) | 0,022 (kp) 0,094 (kp) 0,092 (kp) 0,093 (kp) 
ll (em) | 0,013(kp) 0,079 (kp) 0,085 (kp) 0,087 (kp) 
12 (cm) | 0,007 (kp) | 0,049 (kp) 0,054 (kp) 0,057 (kp) 
13 (cm) | 0,002 (kp) , 0,027 (kp) 0,032 (kp) 0,033 (kp) 
Im zweiten Falle geht man von der Gleichgewichtsbeziehung 
oM, OLen 
a a 
oy 0X Qu: 2) 
(Q,z = Querkraft am Schnitt y = const) 
aus und errechnet für einen Schnitt y = const 
P, P, Rz 
GR oT,, 2 
a ee ee ee Vene 
0 R, Po 
und weiter 
4 2 = 22 
er (Myto) + {Tay)r, — (Tayr,} = EL ge (22). 
Die Integration in y-Richtung erbringt schließlich für ein endliches Plattenelement: 
Y pP 
(M, tot)y ag (M, tot)y=o -r J (Tr, Zur (Tz)r,} dy == CH u (23). 


0 


Am Rande y = 0 ist wegen der Randbedingungen (M, = 0) auch (Mytoı) gleich Null, so daß 
sich vorstehende Beziehung zu 


Y 
j m nz: « 
(M ytot)y + | {(Tıyr, — (Try)p,} dy = ee er (24) 
1) 


vereinfacht. In Bild 17 stellen die eingezeichneten Kreise die Integrationswerte?) für die ver- 
schiedenen Schnitte y = const dar, und die Vergleichsgerade ist entsprechend y- P/2(P = 1(kp)) 


+) Da sich die Spannungen jeweils auf ein bestimmtes Schnittufer beziehen und sich am endlichen Platten- 
element entgegengesetzte Schnittufer ergeben, wurden ebenfalls die (7zy)r,- und (Tzy)r,-Werte mit ent- 
gegengesetztem Vorzeichen in die Rechnung eingesetzt. Der zweite Summand in Gl. (24) beträgt etwa 10—15% 
des Wertes von (My tot)y- 
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eingetragen. Hier ergeben sich z. T. größere Abweichungen, die auf gewisse systematische Fehler 
ne er: eine a Fehlerfortpflanzung schließen lassen. Zu den gröberen Fehlern dürfte 
sicher in starkem Maße der Verlauf der Isoklinen verantwortlich sein, da er in ‚Randnähe nicht 
genau anzugeben ist. Aus diesem Grunde sind die genauen Lagen der singulären Punkte ‚am 
Rande S, bis S, schlecht zu ermitteln. Sicherer sind S, und $,, da bei entsprechenden Gleich- 
gewichtsbetrachtungen in den Eckpunkten nach G. RUMPEL [29] auch das Torsionsmoment Tag 
verschwindet. Relativ sicher sind S, und S, (s. 
Isoklinenbilder). Auf Grund der Randbedingungen 
an den freigestützten Rändern, wonach die Haupt- 
momentenrichtungen unter + 45° zum -Rand ver- 
laufen, und auf Grund spezieller Betrachtungen 
des Richtungsfeldes der Hauptmomente muß zwi- 
schen S, und S, ein weiterer singulärer Punkt liegen, 
der aber spannungsoptisch nicht sicher zu ermitteln 
war. Die Hauptmomente an den freien Rändern 
liefen bei unseren Aufnahmen parallel zu diesen 
Rändern, so daß im Gegensatz zu anderen Ab- 
handlungen keine Torsionsmomente (T,) an den 
freien Rändern auftreten, wie auch statisch ein- 


; y 
Mu N Oedn -(Ty)n, dy (Kp em) 


0 n 2 3 4 5 6/cm)  ]Jeuchtend ist. 
Bild 17. Kontrolle der Ergebnisse nach der sich an einem end- Bei Betrachtung der Momenten-Kurven hat 
lichen Plattenstreifen ergebenden Gleichgewichtsbeziehung sich ein Sachverhalt ergeben, aufden schon H. Vor 


aufmerksam gemacht hat. Danach treten nahe der 
stumpfen Ecken Biegungsmomente auf, deren Vorzeichen entgegengesetzt zu denjenigen in Last- 
nähe und deren Werte relativ groß sind. Das Verhältnis der Hauptbiegungsmomente M, vom 
Lastpunkt und dem Spitzenwert nahe der stumpfen Ecke dürfte etwa 2,5:1 betragen?). 

Die Übertragung der am Modell erhaltenen Ergebnisse auf das natürliche Bauelement 
(Hauptausführung) ist leicht durchzuführen, wenn die Poıssoxschen Konstanten der beiden 
Schichten der Versuchsplatte und des natürlichen Bauelementes gleich sind. Setzt man weiter- 
hin voraus, daß die beiden Platten in ihrer geometrischen Form, Belastung, Deformationen und 
Randbedingungen ähnlich sind, dann spielt nur der Kräftemaßstab bei der Übertragung eine 
Rolle. Da die hier dargestellten Ergebnisse für eine Last von 1 (kp) gültig sind, hat man die 
am Modell erhaltenen Momente lediglich mit der gewählten Einzelkraft P zu multiplizieren, um 
diese Plattengrößen für die Hauptausführung zu erhalten. Dabei muß man natürlich bei der 
Zuordnung der einzelnen Plattenpunkte (entsprechend der Ähnlichkeit in der geometrischen 
Form) den Längenmaßstab berücksichtigen. 

Die Verfasser danken dem Direktor des Instituts für Angewandte Mathematik und Me- 
chanik, Herrn Prof. Dr. K. SCHRÖDER, für die stete Förderung ihrer Forschungsarbeiten. 
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. The Exact Theory er Flat Slabs 


By Kazuvosı Oxo 


‚In der vorliegenden Arbeit wird die Biegung einer gleichmäßig belasteten und durch Reihen von runden 
Säulen gestützten Pilzdecke behandelt. Anfangs werden die Biegemomente in der unendlich ausgedehnten 
und im Ursprung belasteten Platte ermittelt. Dann kann man die Biegemomente in der Pilzdecke berechnen, 
und zwar durch Summieren der von der individuellen Last und von dem individuellen Auflagerdruck hervor- 
gebrachten Biegemomente in der umendlich ausgedehnten Platte. Zu diesen Berechnungen sind die Weier- 
straßschen elliptischen Funktionen benutzbar. 


. The present paper deals with the bending of a flat slab which carries a uniformly distributed. load over 
its entire surface and is supported by rows of equidistant columns with circular sections. In the first place 
the bending moments in an infinitely extending plate which carries a load at the origin are solved. Then the 
bending moments in the flat slab are obtained as the sum of bending moments in the infinitely extending plate 


caused by the individual load and by the individual reaction from each column. The Weierstrassian . 


elliptie functions are used conveniently for this purpose. 


B nacronmeii pa6oTe paccMmaTpuBaetca IIPoruO PABHOMEPHO HATPy;KEHHOTO ÖC30AJIOYHOFO 
epeKpbITUA IHONIMEPTOTO PANAMH KPyTiIbIX KOJIOHH. CHayana ÖyAyT HOJIy4eHbI U3TUOAIIIME 
MOMEHTBI ÖeCKOHEYHOÄ ILIACTHHKH C HATPY3KOÜ B Hayallıe KoopnmHaT. Ilocıe 3TOTO MOHHO 
BbIYUCHHTB MU3THÖAIIIMEe MOMEHTEI B 6E306AJIOYHOM HEPeKPBITUM, CYMMMUPYA M3TNÖAIOINMe MO- 
MEHTbI B ÜeCKOHEYHOH ILIACTHHLbIE, BBIBBAHHOE COÖCTBEHHOH HATPY3KOlÜ HU COÖCTBEHHBIM ONOP- 
HBIM MABJIeHMeM. lH 3THX BbIyucsIeHNÜ MOSKHO HMCNOAB3OBATB IANIUNTMYeCcKHe PYHRUMH 
Beüepmrpacca. 


Prefaee 


We know two ways of solving the problem of a flat slab which is subjected to a load distri- 
buted over its entire surface. One is represented by Lewe’s solution!), by which the deflection 
and the bending moments in the slab are represented by double trigonometric series. This 
solution has the defect that the series which express the bending moments converge too slowly. 


The other solution was developed by S. WoINOWSKY-KRIEGER?). The slab is divided into 
a number of support-strips (Stützenstreifen) and field-strips (Feldstreifen), and the deflection 
of each of them is given by a single trigonometric series. The coefficients of the terms of these 
series are obtained by troublesome calculations from the boundary conditions of these strips. But 
this solution has the advantage that the series which express the bending moments are rapidly 
convergent. 

We have another solution?), which, ignoring the dimensions of the columns, gives the bend- 
ing moments in the slab with ease. But the resulting errors are too large to be disregarded in 
the slab close to the supporting column. 

In these solutions and in other solutions which have been hitherto proposed, the deflections 
of the slab are always expressed by single or double trigonometric series, but it seems advisable 
to use some other functions more suitable to represent the deflection and the bending moments 
in the flat slab. We can easily solve the deflections and the bending moments in circular plates 
carrying symmetrically distributed load by determining the constants of integration in the 
general solution of the fundamental differential equation®) expressed in the polar coordinates 
so as to satisfy the boundary conditions. In this paper we shall denote by the symmetrically 
distributed load the distributed load with intensity, which is expressed as the function of the 
distance from the center or origin, but not influenced by the direction from the origin. 


If we assume that the flat slab is a circular plate of infinite radius carrying a uniformly 
distributed load and receiving equal reactions from rows of equidistant columns, we shall be 
able to calculate the bending moment in the flat slab by a process similar to that for the problem 
of the circular plates of finite radii. The author has developed this idea and has represented the 
bending moment acting on any section of the flat slab by a single series which converges very 
rapidly. This solution can also be applied to the flat slabs supported by columns arranged at 
the vertices of parallelograms. 

1) Lew, Die Lösung des Pilzdeckenproblems durch Fourizgsche Reihen, Der Bauingenieur 1 (1920). 

2) S. WoINOWSKY-KRIEGER, Beitrag zur Theorie der Pilzdecken, ZAMM 14 (1934). 

5) A. NAvar, Über die Biegung durchlaufender Platten und rechteckiger Platten mit freien Rändern, 
ZAMM 2 (1922). 

4) A. NApar, Elastische Platten, 8.55, Gl. (13). 
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Chapter I 
Bending of an Infinitely Extending Plate under the Action of Lateral Loads and Moments 


1. Bending of a Plate under a Symmetrically Distributed Load 
We shall begin the theory of bending of plates using the following notations: 


za : Rectangular coordinates. 
r,6 : Polar coordinates. 

20, : Oblique coordinates. 

E : Modulus of elasticity. 


v : Poisson’s ratio. 

h : Thickness of the plate. 

D : Flexural rigidity of the plate. 

p : Intensity of the distributed load. 

w : Deflection of the neutral surface of the plate in the downward direction. 

M,M, : Bending moments per unit length of sections of the plate perpendieular to the x- 
and y-axes respectively. 

My,» My.: Twisting moments per unit length of sections of the plate perpendicular to the x- 
and y-axes respectively. 

s : Radius of a circle in which a distributed load acts on the plate. 

2a,2a, : Distances between the columns in the u,- and u,-directions respectively. 

: Total magnitude of a distributed load acting on the plate inside a portion bounded 

by a circle of radius s or bounded by a parallelogram with sides 2 a, and 2 a,. 

'% : Magnitude of a moment acting on the plate. 


Using these notations the deflection of the plate can be determined by integrating the 
differential equation: 


02 02 2 ) | 
(+ )9=5 an ee Sr 
or | 
PR ER, 
6 Eu a .o an) w— Tr TOM BMI SEPEaEL BE (1.2), | 
where 
ke Eh? 
> EI) I 


When we use the recetangular coordinates, the bending moments and the twisting moments 
are calculated by the equations: 


0? 2 
age 2 +» = 
Or Oy® 
e?ıw o2w 
MD a Een 
2. ni erw 
M;y =—M.=(l—»)D dr öy 


Ihese moments have the positive sign when they are directed as shown in Fig.1. 


Introducting new notations 


D /o:w o:w 
A JLcEn 
ur D o:w 7) 
M,=-—5 (nl RU N. 
Myx Qy dx M, En Rn) 


Fig. 1 


7 dedy 
M, M, and M,, can be represented in the following forms: 
M, = (1 -- v) M, 1 (1 — >») M,, 


M‚=(l1+)M—(-»)M,, 0 ER ne EEE 
My=—M,„=—(l —»)M, | 


> 
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—: 


The bending moments M,, M, and the twisting moment M,;, are determined easily from 
M,, M, and M,, so in the following part we shall calculate M,, M, and M, instead of M,, M, and 
M;y. When the deflection w is represented by the polar coordinate, M,, M, and M, are cal- 
culated by the following equations: 


M--5[% 1 ow a) 
2 \or r or r2 962 
nl : en Ge) es 20 
ar E EL = R = sin 2 0, We (1.7). 
me fr nt ? 


1 &w 1 ow 
r or 00 r? 00 


+2[ )eos2 0] 


Now we shall discuss the bending of an infinitely extending plate when a symmetrically 
distributed load is applied to the plate on a portion bounded by a circle of radius s with its center 
at the origin. The distributed load is assumed to be expressed by the following form: 


Tr 2 
irarle) a ne | BL 5 (1.8); 
for r>Ss | 


where A is any arbitrary real number, not necessarily integer, but satisfying the condition: 
1>—2. 


The total magnitude of the distributed load becomes 


P=2n2 3 75 ‚Hat ATS (1.9). 
A 


Substituting equation (1.8) for p in equation (1.2) and solving this differential equation, 
we have 


w= =. je: +. 8) log) tn: | for es, el 
/2) S A E n Ss 

u-..p|(! & + log @ t Go —% 8? +0, 5? log 7) 

si Pı ip Nr ; 
+ D 2 0 £ 9)2 (2 = 4) = for I = A Man (1.14), 

and in particular for r — 0 
BD t 2 2 e Ss c 

=D (0a + a, s?log 7) DET (1212): 


where C, and l are constants of integration and x, &% and a, can be calculated by the expressions: 


bern > ee » 


Pi 


} 
5 
»P=2ns 2, 4,9% TRETEN AR) 


EN WB 
Kz Pi Dr S= > 0 4% 
4 


When the radius of the plate is finite, the constants C, and ! can be determined by the 
boundary conditions; but as we have no information about the boundary conditions at r — © 
for the infinitely extending plate, we shall leave the constants indeterminate. 
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The bending moments are obtained by substituting equations (1.10) and (1.11) for w in 
equation (1,7): 


M=—;,(1t18 r)- 
r % 5? 
er (' y; se a a (1.14); 
P &$ 
Mm ;,[1-) in 29, 
for res 
P ru s2 P; er TER 
M=-—,, (12 + log 1) ge I 
Be 5? pi au 9 
MT er Tea 9 GE re 
aa 
re Peer. ) DR 
for Tu 
and in particular for r = 0 | 
DB s\ { 
| | Bel. 


M=%, M,=% 


Values of &,, &, and «, for various modes of distribution of load of which P = s?, are given 
in table 1. 


2. Bending of an Infinitely Extending Plate Produced by a Moment with its 
Axis Parallel to the Plate 


We shall discuss the bending of an infinitely extending plate subjected to a moment, re- 
presented by a distributed load applied to a portion of the plate bounded by a circle of radius 
s having the center at the origin: 


A—1 
p= 24 (*) sin (d — y) for tS3, 


PpP=0 for r>s 


. (2.1), 


where A> —.2. 

The total magnitude of the distributed load is of course zero, but 
the load produces a moment with its axis directed to 0 =y (Fig. 2) and 
with the magnitude 


ng. 2 q | 1 DR 
“ True Die ee 


Substituting equation (2.1) for p in equation (1.2) and solving this differential equation 
we have 


a Sn a 
w— ) \ (' H 2log ‚) + 37 | sin (d — y) for DES ar oe 
ze, 4 | RB ad, Ps N 
v= g7D I — 2, + 2log ı) r-+ ys | sin (0 — y) 
st (do {F 8 r y PM se 
75 16 (5) log Ak; 2 Ia+MÜH 55 ) [sin (0 —y) “2 
for ER 


where 2” denotes the sum of the terms in which A # 0. We can easily show that the value of 
A 


the constant of integration ! in equation (2.3) must have the same value as lin equation (1.10), 
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Table 1. Values of &,, x, and &3 


Tue - nA * 
> - - „ J r 3 art; En u 
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provided that the boundary conditions at r > ©© are equal to each other. The constants ß}, Pa 
and /, in equations (2.3) and (2.4) are calculated by the expressions: 


DD re 
AT=as I are nn (25). 
armer) 


By substituting equations (2.3) and (2.4) for w in equation (1.7), we have the bending 
moments 


a) h 
er 
m cos (d—y)sin20 As 2,6 
M, = ee v 23 ae er 
T muznm20 Br en 
M = In | wre 2 r 


for TENS®, 


M, =2M,sin (# —y), 
M, = — M,sin (0 +y) — M, sin 30 —y), ee 
M, = M,cos (9 +%y) + M,cos (39 —y) I 


for =S, 
a: 1 r ‚q 2 Aryı 
% de u a I 72a | 


str q ri 
M=-7 3 ar9arS\s) 


The values of f,, ß; and ß, for various modes of distribution of load of which T=xs° 
are given in table 2. 


where 


. (2.8). 


The diagrams in Table 2. show the distributions of load along the radii with directions 
7 7 
Per, and (Fiers- 

We have hitherto assumed that the moment 7 is acting on the plate in the form of distri- 
buted load expressed by equation (2.1). But when the moment 7 is applied to the plate from a 
column of radius s supporting the plate, we should rather assume that the ceircular portion of 
the plate of radius s immediately above the column is absolutely rigid. It is also assumed that 
the edge of the plate is clamped along the circle of radius s. Hence the following boundary 
condition must be satisfied: 


ow w 


Pk for r=S$S I 


Substituting equation (2.3) in this condition, we have 


Fe Re... 


3. Deflection of Plates Carrying Several Loads of Various Kinds 


We can obtain the deflection of a plate carrying several loads of various kinds by super- 
posing the deflection due to the individual load. Equations (1.10) and (2.3) show that the deflec- 
tion of the plate increases infinitely for r— 00 but we can make it to approach zero for r— oo 
by superposing several loads of which the kinds and the magnitudes are properly chosen. In 
this article we shall show several examples. 


re | Table 2. Values off, Ba and ß, 


Pl 
p=lm (4-+2) =) sın (d —y) 
i>—2 6 


B, = 0 


(= 


=: Bs = —3 


p=4— sin d—y) 


2 
Beth 


Ji—1 
p=iim ara [7] nen) 


A>an 


=! R=0; 
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(i) A moment and two loads on a plate 
When a moment T acts on an infinitely extending plate at the origin with its axis 6 A 
we can make the deflection to approach zero for r — oo by supplementing the moment with two 


= 
loads of magnitudes P and — P acting at the points [ =a,09=y— n and [r =4,9=y+ ) 


respectively (see Fig. 3). ‚ ] 
The deflection of the plate is obtained by superposing 


pP equations (1.10) and (2.3) as follows: 
T r Pi: 
+ Bi (ri + 2) log! +6,+%% 
BED em Bars 
P (24 Qlg?+C, + GN) 
—EnD Leere ie 
where 
R=r+2rasin0—y)+@ | 
R=r—2rasn0—y)+@ | 
From equation (3.1) we see that w approaches zero for 
r — 00 regardless of the value of 8, when the following condi- 
tion is satisfied 


aa 


T=3u Pi. 20,02 En 


(ii) Four loads 
When four or more loads act on an infinitely extending plate at points situated on a straight 
line on the plate (Fig. 4), we can make the deflection of the plate to approach zero for r— oo 
regardless ot the value of by selecting the magnitudes of the loads and their points of application 
so as to satisfy the conditions: 


B R 
5 B mn ru =Uu >P,5=-0 8% 
F; KK a Er wiere al. 
2 In these equations x, is the abscissa of the point 


of application of the load P„ on the assumption that 
the x-axis of the rectangular coordinates is drawn so as to pass through the points of application 
of the loads. The values of a, s? should be equal for all the loads when they are represented 
by equation (1.8). 


(ii) Three moments 


When an infinitely extending plate is subjected to three or more moments with their axes 
parallel to one another and acting at points situated on a straight line on the plate which interseets 
the axes of the moments at right angles (Fig. 5), we can make the deflection of the plate to 

approach zero for r— oo regardless of the value of 9 by 


0 x x, x ; re : 3 
5 ] 2 3 selecting the magnitude of the moment and the poin 


- - = application so as to satisfy the conditions: 
1 1: J ı m 1. 
r-neß; ZT, um= 0,7 BE 
where n=1, 2...» 
(iv) Many loads and moments 


As a general case we shall find the conditions which are used to make the deflection to 
approach zero for r > oo, when arbitrary number of loads and moments are acting on an infinitely 
extending plate. We shall denote the magnitudes of the load and the moment which are acting 
at point (X, Y„) by P„ and T,„ respectively (Fig. 6). Then we have 

EPr=0, 

& PR, — 2 T,sinyn—=0, 

2 PaYyn + 2 T’7„ C0S Yn =, 

2 PR 22 —22 T,siny=0, | BT 
2 Pnyn +2 2 Tn Yn cos m =0, 

2 Pain Br (X, COS Ya — Ya Sin yn) = 0 


4 
I 
W 
b 
j 
I 
\ 
BE 


I 
S- 
1; 


=; 
&. 


\ ERS ER f R NWaRII 


ERNNERE 


om 


DEN > 


\ 


EN y — re 


: Mn 
BR. YA } f m 
ER & \ f 
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Generally speaking, when many loads or moments are acting on an infinitely extending 
plate, we can divide them into groups, each consisting of such loads and moments as satisfy the 
condition that the deflection of the plate at r— oo approaches zero, in other word, it is O(1/r°), 
where r is a distance from the group and e is a positive number. - 
When infinitely many groups of loads are acting on the whole sur- 
face of the plate, we must note that the superposed deflection of the 


plate does not always converge even if the deflection due to each (Xn » Yn) 
group is O(1/r*). But in the present case we can easily verify that BNA 
the sum of the bending moments produced by these groups of loads ; 
converges, because the bending moments produced by each group > 
become O(1/r?*®), as a result of differentiation of the deflection with Me 
regard to vory. Fig. 6 


4. Bending of an Infinitely Extending Plate under Loads and Moments Acting on 
a Portion of the Plate Bounded by a Parallelogram 


We shall consider the bending of an infinitely extending plate produced by uniformly 
distributed load p applied to the plate on a portion bounded by a parallelogram with its center 
at the origin. As shown in Fig. 7 we shall denote the lengths of the sides of the parallelogram 
ABCDby 2a, and 2a, and the angles between the x-axis and the sides of the parallelogram 
by y, and y,, provided that 


Div Pi en ui Bar). 


Using oblique coordinates u,, u, having their 
directions parallel to the sides of the parallelogram, 
we obtain the deflection of the plate at the point 
(R, 6) by integration 


w= ul j (g.108 7 + c) du, du; sin ys, 
Ber 4.2), 
where g2 
e= Rr—2Ru,cos(d—y,)—2 Ru, cos (y;— 0) 
+ +3 +2u,U,C0Sy, 
Fig. 7 
A Pe ee Pr BE Fr (4.3). 
For the point which satisfies the conditions 
R>a, I (4.4), 
we have 
P a? + a: R a+a 
ir RR = 2 — +6+--— 
e BD ılR au De 108 es 
a Ü (8 6052 9 +0[(- 4.5 
+7 c0s2( 7 Male 208 (Y—0) + ve BESTE ee (4.5), 
where 
ER SR A DE Ve a rn (4.6). 


When the reaction P acts upwards over the inner portion of the plate bounded by a circle 
of radius s with its center at the origin, we obtain the deflection of the plate by changing the 
sign of equation (1.10) 


P 


R 
en (R? + a, 82) log T RE (4.7). 


The distributed load in the parallelogram and the reaction at its center are in a state of 
equilibrium, but the sum of the deflections expressed by equations (4.5) and (4.7) does not 
approach zero for R— oo. For the purpose of letting it approach zero we shall add moments Ir 
and T, per unit length of the sides of the parallelogram with their axes parallel to the sides BC, 
DA and AB, CD respectively and two upward forces with magnitude Q at vertices A and C 
and two downward forces with magnitude Q at vertices B and D (see Fig. 8). The deflection of 
the plate at point (R, 6) produced by these moments and forces is expressed by the following 
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equation, provided that expression (4.4) is satisfied: 


mi (7,4 Tysiny + Q sy) [1 41087) 


w=—. 


m 1 
(nnd) + z0 0 nn) + (| . (48). 


From the condition that the sum of the deflections expressed by equations (4.5), (4.7) and 
(4.8) approaches zero for R— oo regardless of the value of 0, we obtain 


= Br 
1, = sinn — a8, e 


2 
= Persien — ns, (49). 


(a8 E) 
0 Ban sin Ya COS Ya 

We shall hereafter denote the above-men- 
tioned superposition of the uniformly distributed 
load over a portion bounded by a parallelogram, 
the reaction acting at its center, the moments 
applied along its sides and the forces acting at 
its vertices by the term “fundamental group of 
loads of the first kind”. 

We shall next consider the bending moments 
at any point (r, 6) on the plate produced by this 
fundamental group of loads. The point (r, 6) may 
be inside or outside the parallelogram in so for as 
it is outside the circle of radius s. 

The bending moments produced by the reaction P are obtained by changing the sign of 
equation (1.14) 


e PR 5 
M, = ye ( =—- r? ) cos 2 v) 3 a Pe A (4 10) 
RP SA . 
= nn 2 
M, = \; )sin20 


We can express the oblique coordinates u,, us, by r and ®, and the polar coordinates r, ® 


by u, and ı,, 
BE a | 
sin Ya 


in Bu are (4.11), 
w=r — 
SIN Ya 
r=W +W +2uU,C0Sy, | 
d tan. SU yı Fun Ein ps (4.12), 
U, COSYy, + U, COSYyg | 
where 
Yy—n<d<y for 10% 
Yy»<l9<y+n for u. 0, 
a er ls leer (4.13) 
I=y—n for w=0,,<0 


Using equation (4.12), we can regard the bending moments expressed by equation (4.10) 
as functions of u, and u, and express them in the forms: 


a y 
M= I {Fi(u,, ug) + 0, kK,l(u,, u,)}, | 
Te 
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where 


7.58 
= Ss Ep PC 2 BE Tel Te LS N 6 Re (4.15), 
Dis ArOSSIN yon a (4.16); 
Fl, u) = 1 +10, 
cos 29 
Tel) a NE Ner Bee (4.17); 
sin 2 
F;(u,, u) = - 2. 
K,u,u)=0, 
S 26 
Ku, u) = — 5-: = » 
EB 7 N (4.18) 
S sin20 
K,(u,, U;) = — IE Te 


In the following part we shall use the abridged notation (+ -+) fe + a, y+ b) for the 
expression 


«+ y+b)—-a«@—ay+b+fa—ay—b)—-f@+a,y—b)... (4.19), 
and as a particular case (+ -+)fz + a + b) for 
e+a+b)—- f«@—ar+b)+fz@—a—b)—-f@+a—b)..... (4.20). 


Then we can express the bending moments produced by the uniformly distributed load 
over the inner portion of the plate bounded by the parallelogram ABCD in the form: 


P 
Me T, Fe) (U, 0, UzseQ,), | KOREA Te (4.21); 
N I | 
where 
sin Ya 
28 


G,(u,, U,) = 


{(uf + u3) cosy, + 2 u, u,} 5 — log =) 


+ {120 + (0 —w)}siny, 


9 


1 
zul (cos2y, + cos 2y,) 


sin ya 
(u, u)= 5 © | 


+ f(- 182 sin2y, + w sin 2y,) log . 


2 aA), 
+ 1409 cos 2y,—.u3 (d — u) cos2y,}sinyz| , 
Gi, U,)-= ne = e: u, u, (sin2y, + sin 2y,) 
+ je cos 29; — u2 cos 2 y,) log T 
+ 26 sin 29, — u? (9 — u) sin 2 Z sin ya 
RN, my reis tn,] (4.23). 


uam yore f ee 


The bending moments produced by the moments T,, T, applied to the plate along the 
sides of the parallelogram ABCD and the forces with magnitude Q applied at vertices A, B, C, D, 
are 


r 2 
M; = pe ++) Ha aW+ a)+ Ne ük+-+) Lu Ft Wa) . (4.24), 
?==1,2,3% 
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where 
sin r & 
Hl u) = "gg (+ eos7alog 7 


+ 0 sin y — a3 (0 — u) sinn); 


a, sin yz [,. Lo . 
H;(u,, U;) = 54 = [sin 2 y; log ; 0cos2 2, 


nn, Qz SIN ya [sin 2y, log T — (9 — u) cos271] 
a 


1 


a 


dy Az 7) ne 
+5 ae (28 -+9)+ 18 a, cos (2 Yı) 
H,(u,u)=— ei [eos 2 y, log T +9sin 2 2 
3 
+ er [eos 2 y, log 7 + (# — u) sin 2 2 
1 
re 20 +y)+ Ban 8 —y;) J 


1 
L,(u,, u;) = — An Ks 


Ls(u,, u) = = [ein 2y, — sin 27y,) log 5 — (d — u) cos2y, + 9 cos2 7! i . (4.26). 


Liu,n)= => | cos2y, + cos 2y,) log - — (9 — u) sin 2y, + ® sin 2 2 


None of the bending moments expressed by equations (4.14), (4.21) and (4.24) approaches 
zero for r— 00, but we can easily show that the sum of the three bending moments M, expressed 
by these equations is O (1/r*), and the same can be said about M, and M,. 

Now we shall regard the surface of the plate as a plane of complex numbers, and any point 
on the plate can be represented by a complex number z where 


j 
i 
t 
| 
H 


z=r(cos0 +isin ®) | 
= U, (c08Yy, + isiny,) + u (cos yy + I sin y,) j 


The lengths and the directions of the sides of the parallelogram are also expressed by 
complex numbers 2 w, and 2 o,, where 


2 ©, = 2a, (c0osy, +isiny), (4.28) 
200: = 2, (eos yo ana JE ur: 
Here we shall introduce new notations & and £,, and let 
eT U, R-J 35) 
ee 3 3() q 
VD; 
(4.29), 


and regard them as functions of z, 


Then the functions F,(u,, u,), G,(u,, u,) and H,(u,, u;) in equations (4.17), (4.22) and (4.25) 
respectively can be expressed as functions of the complex variable z: 


Fu, u)= RE, 2 (4.30) 
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where 
Be ee ee EN St: 
El N ee (4.32), 
where 
ee 
= + zn : ee N (4.33), 
22 1 zZ 
G(z) = 8010, > = log 7) Ba ae ge (4.34); 
Ku. 0) SRH ea (4.35), 
where 
Bea (4.36) 
2 ee .36), 
ale, 05 2 
H() = A m %) N ee (4.37). 


In these expressions and hereafter $ (log z/l) is defined so as to take the following values 
provided lisa positive real number: 

#4 

| arg — and (108 = Te ee ie (4.38). 


The value u in equations (4.33) and (4.36) are defined to take the following values: 


w=0R or arg ma <% (log 7) < arg © +2, | 
e | a (4.39). 
= 200. for ago, Ha <A[log 7) <argn +” | 


Now we shall consider the relations between the functions F,, F,; G,, @3; H,, H, and the 
functions F,, G, H, respectively. Using &, and £, in equation (4.29), we can express z in the form 


DEE OR ee (4.40). 


Substituting this expression of z in equations (4.31), (4.33) and (4.36), we can represent 
them as functions of the four variables £,, &,, &, and w,, which are expressed in the forms 
Fu) = Fol&ı, &; © @;]1 
a (4.41). 
Ho) = Holdı, 835 @1, @5] 


We can easily show the following relations to exist between equations (4.17), (4.22), (4.25) 
and (4.41): 


A F(u, 4) = R{Fyl[E 5 @y @]}; ) 
, 1042,=0 ER 2 ,0(4.42); 
— Fz(u, u;) + I F,(u, U) = — 5 G don 4 On) Fol&ı &5 @1 @s] | ar 
G,(u,, U) = RGolEı, &; @1 @s]} » | 
. 1 IR) rd 4.43); 
— Gi 0) +1 Ol 1) = — 5 (Br za; + za) Oli no | 
H,(u, 4) = R{Hol&, 8; © ®s]} ; | 
ee (4.44); 
— H;(u,, u;) + i Hu, u) = — 5 G dar + Da) Hol &35 @1, @;] j 
where ®, and @, are conjugate complex numbers of ®, and w, respectively and are 
a ET STR (4.45). 
OD; = A, (C0S y; — sin y;) 
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Now we shall express the functions K,, K, defined by equation (4.18) and L,L, defined 
by equation (4.26) as functions of z: 


— K;(u,, u,) + i K,(u,, u,) = K,(2) es re Re a (4.46), 

where R 
K,(z) = ek ta an u ie IRRE ee (4.47); 
— Lf(u,1) till, 1) =Ll) -» -» 0. (4.48), 

where = 
© nn 
Lo) =1L)— Pen = A (4,.49), 
. 

Dr oe a Re 4.50). 
L(z) a; Sn o Pas I ( ) 

Chapter II 


Bending of a Flat Slab under Uniformly Distributed Load 


5. Representation of the Bending Moments in the Form of Elliptic Functions 


In the preceding article we have considered the bending moments produced in an infinitely 
extending plate under the action of a fundamental group of loads of the first kind. If we divide 
the plate into infinitely many parallelograms by drawing two sets of equidistant parallel lines, 
and assume that the same fundamental group of loads is 
applied in each parallelogram formed by the lines. 


Then the plate represents an infinitely extending flat 
slab supported by columns arranged at the centers of the 
parallelograms and carrying a uniformly distributed load 
over its whole surface (Fig. 9). The moments and the forces 
which are assumed to be applied to the plate along the 
boundary of each parallelogram and at its vertices re- 
spectively, are counterbalanced by those acting along the 
boundary of the neighbouring parallelograms, and we can 
assume that no moment or force is applied on the slab flat 
except at infinity. 

By the use of the oblique coordinates, the positions 
of the columns supporting the flat slab are expressed by 
(2m a, 2na,), where m and n denote positive or negative 
integers or zero. The bending moments in the flat slab at any point can be obtained as the sum 
of the infinitely many bending moments, each of which is produced by the fundamental group 
of loads acting on each parallelogram. As the bending moments expressed by equations (4.14), 
(4.21) and (4.24) consist of two terms, one of which is independent of a, k and the other has the 
coefficient a, k, the bending moments in the flat slab can also be expressed as the sum of two terms ‘ 


M;=M+Yo,kMj, Tea Dee 


where 


P 
ER N ‚ - 
4 = {Fu —2ma, w—2na,) 


mn 


+=+: +J)%w—2ma + a, w—2naH+ a, +0) 
+: +) Hy —2ma+a, w—2na+ Q;)} 
Ei 
M; = yir= r {Ku —2ma, w—2na,) | (5.3) 
++: +)Liw—2ma +4, w—2n A; + Q,)} | 


for 0 2 s, where 
= (1 — 2m a) + (yw—2na,) +2 (m —2m a)(w—2na)cosy .. (5.4) 
As mentioned in the preceding article each term of the series in equations (9.2) and (5.3) 


: . 1 
approaches zero with the order of — , when o approaches infinity, we can easily show that the 


bending moments represented by equations (5.2) and (5.3) are convergent. 
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We see here that the bending moments M,, M, and M, are expressed by periodic functions 
of u, and u, with periods 2 a, and 2 a, respectively. 


Substituting equations (4.30), (4.32) and (4.35) for F,,G, and H, respectively in equation 
(5.2), we obtain 


P 
ee 2.6.5) 


for 


where 


- en 
ed =ailg—htäh)t ZFle— Om) Ar mal —m +) ee 


I GEN nF 020) air Ee) IE ne ER) 


and 
DEN 20T ee ee (5.7): 


The first term of the right-hand side of equation (5.6) and the second term in the double 
braces are added to this equation in order to make the imaginary part a convergent series and 
to make the function exp {D,(z)} a continuous, odd and modular function for all values of z 
[See equations (6.1) and (6.33)]. If we denote by m, and n, integers which satisfy the conditions 


1 1 
A757 <m SH ee 
1 a (3.8), 
air Dr 
Ön,n and e in equation (5.6) are defined by expressions: 
Oh loss Em m ans 
en, tor mMem, Nenn) = 
On N RER ER a ca 
a) oe 
1 
Bo OL em, | 
ee ae ee ee (5.10) 
ee; for sam | 
Substituting equations (4.33) and (4.36) for G, and H, in equation (5.6) respectively, we 
obtain 
D,c@) = De) +ri ae u ae re re. (9-18); 
a 0] 
where 
Bl 2 , g\ @z 
Bl) = ni or +ri G +n ni) Ar 
Br 2 — An 
+ 2, [Fol — Amm) — I nn me e) RIM 
hack ) Me I mm + © +) 
Ed  Hagia | 
Now we shall discuss the function ©(z). Differentiating ®(z) three times with respect to z, 
we have 


DH) = I a Im) + EIER Im 39H @) | 6.13) 
+(+:+)H"@— 2,4 014 a9} 


‘ 
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where 9 
Ro)=-3 ] 
—ı 
en, 5.14). 
O-TTRE - 6.14) 


DEN Le 1 © @g 
Ne ak 7 = 


Substituting these expressions in equation (5.13), we have a 


vi 
Dr Das 


x sh. BEE ee 
ud ee En a 


1/0, , 0%, 2 
er wle. +) ER e J 
This equation expresses the sum of the groups, each consisting of several terms. As each 
term approaches zero with the order of (€ — 2,,,)? when |z— Q,,„| approaches infinity, we 
can represent ®’’(z) as a sum of several convergent series. It is easy to show that these series 
are doubly periodie functions with periods 2 », and 2 w,, so we can express them by WEIER- 
STRASSIAN elliptic functions 


= — hd. 


4 0, ®z 


(+-+J)!@+o, + ®,) 


1 To; as F 
unse -+)P@ + + ©) 
‚Remembering the quasi-periodieity of {(z) and the periodicity of $’(z), we can show that 
the right-hand side of equation (5.16) reduces to zero except $’(z) and obtain 


Dipl 0 m (5.17). 


Integrating the above equation three times, we have 


Nuke o(z) 2 
®(z) == log Im, En o(z) ee N (5.18), 
a p(z) is a quadratic function of z. Substituting the above equation in equation (5.11), 
we have 
z 


20, 


De) lg 1 e) +trie, 


20, 


‚  Fromthe condition that M}, or P/4r- R{®,(z)} hasthe periods 2 o, and 2», we can deter- 
mine the function g(z) 


Ko) „2 
2 


A (5.20), 


where C, is a constant. Substituting equation (5.20) for p(z) in equation (5.19), we get 


erg 
(2) = log 3, Fe +ri 3 Zw (5.21). 
The value of C, can be determined by caleulatin i 
R g the values of ®,(z) for a particular value 
ae seen ni 1) ac En and then by equating them. Equation (5.5) Reis not be applied 
‚alue ol z where |z) < s, but we shall assume that ©,(z) is also defined i 
for the value of z where z approaches zero, then we have ea; a 


lim !®,(z) — log — = log2r + 210g Q,(r) + = Or (5.22) 
20 2 | 6 Br. 
where 


rund rn ER UN 2 (er 


Pe N 
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and 
Q,G) = II“ — exp (2nzir), | 
or E ER E RTEL (5.24). 
Ed Pe TasppnreveaT | 


From equation (5.21) we get 
z 
lim [90 — log 3 — (Var Sr (5.25). 


As the values expressed by equations (5.22) and (5.25) must coincide with each other, 
we find 


C, = log2r + 21log Q,l) + = ee (5.26). 
Substituting equation (5.26) in equation (5.21) and representing o-function by 9-function, 


we have 
rn 


D,«) = log d(u|r) — log EU BE, DER Dr MAN (5.27), 


-- where 


: 9, (vr) = 2exp(zwir)sinao—2esp(rir)sin 320 +2exp (Frir)sin INU-— ++» (6.28), 


4 
z 

u ROT ENT (3.29), 
og ei Es (5.30). 

Substituting equation en in equation (5.5), we get 

0 
Mr = — a log 9, (v|r) un 2 2 Led ee (5.31). 

Substituting 
eg ad ER (5.32) 
© 


for v in equation (5.27) and &,/o, for T, we can represent D,(z) as a function of &,, &,, ®, and @,, 
and we shall denote this function by ©,[&,,&; ©, ©z] 


©; | | DO; 


1 
Dulin&; © ©s] = 10g 9, (& +& = + 3 log 27 


on . 6.33). 
— a log ö 0 2 +2 is(& +5 —] | 


We shall next proceed to obtain the bending moments M% and M%. Substituting equations 
(4.42), (4.43) and (4.44) in equation (5.2), we have 


are 
— Mi +iMi=— S (50; + 9.) |Fols—m Sn mir] 
Fr 21 1 
1 1 Fi 
+4 Ha mtz sont zin . 6.34). 
1 1 i 
+(+- +) am +, ont 5: on] 


Remembering that the first term of equation (5.6) and the second term in the double braces 
become zero by differentiating with respect to ®, or @,, we obtain 


u Po EN ER (5.35). 


BZ 
— M* +iM!=— u Terre 0,5, 
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For any function of variables v, 7 we can show 


-;[@ 2 +5, 2) fv, ı) = r Pr le Et >)1@ De | (5.36). 


Therefore we have from ae (5.27) and (5.35) 
A r 1 i i 
— Mi +iMf=7- . a ale + +5 .) los Bol) — z 108 d0) +a ik, 2 (5.37). 
For numerical calculation we shall express equations (5.31) and (5.37) in the form of trigo- 
nometric and hyperbolic functions: 


Pt o 2 
log (2 sin) Hair. +8) —2 >- cos2nzv & (5.38), 


M= ER „= niexp (— 2nzir) — 


—M: HiM=In Da is c00n0—(z +4) 
sin2nzv 
+4i8, PR er eye (5.39). 


© cos2nrzv 
Pre een, 


The series in these expressions converge, when the absolute value of £, is less than unity, 
and as £, approaches zero they converge more rapidly. Generally speaking sufficient accuracy 
is obtained by taking only several terms of the series when the absolute value of £, is equal to 


1 
or less than DE 


As M}, M%$ and M% are periodic functions of z with periods 2 ©, and 2 ®,, the bending 
moments at any point on the flat slab can be calculated by the use of these equations. 

Now we are in a position to calculate the bending moments M/ expressed by equation (5.3). 
Substituting expressions (4.18) and (4.26) in equation (5.3) and remembering equations (4.46) 
and (4.48), we have 


nn 


M, = Am PREISE (5.40), 
’ . ! P - 
—M; +tiM, = FE ED Re ME e (5.41), h 
where 
Pd) = I {Re — Om) + + HL An tatw)} - - - - 62). 


Mn 
Substituting expression (4.49) in equation (5.42), we have 
0) 
FIN en ” 1 
v; (2) a P(z) Fr, 05} 
20 
where 


Y(z) So 3, {Ko er An) al ea £ Je) L(z ro Don == © Ex @3)} ir ee - (9.44). 


mn 


Differentiating this equation, we get 


’@)= = 2 {Ko 02 — Dmım) + es ) L’(z Pr N A + en + ,)} Re: (5.45), 
where 
CR S 
K,(@) a pe 23 ’ 
“ . (5.46). 
Lie) 


I 0, wm 2 
Sııt n r a "aacı ji 3 1 at 
Substituting expression (5.46) in equation (5.45), we have 


/ 5 1 ee, Le) 
w (z) = — > Fe ig Mn Er 
It Mn (2 — N, + (z — 2° ee — BE (w? Eu oe) (z er 0 ” 74 (0 — od)? (5.47). 
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R en equation can also be expressed by WEIERSTRAssian elliptic functions with periods 2 @, 
and 2 @; 


Fe) =5,- + — (+. s+ER+m + a) een (5.48). 


ren 


Remembering the quasi-periodicity of £-function, we see that the second term of the right- 
‚hand side of equation (5.48) reduces to zero, and we have 


, S 0, 
DR = 5: N ee Kap BER (5.49). 
Integrating this equation we get 
S 
Pd a Dede ne A ee (5.50), 
where C, is a constant. Substituting this equation in (5.43), we have 

AS © 

RL) = In 9(z) => Io, + (3 ne En oe (5.51). 


Equation (5.41) cannot be applied for the value of z where |z| < s, but we shall assume 
that Y,(z) is also defined by equation (5.42) for the value of z, where z approaches zero, then we 
have 


1/o. m} ©. 
1 EL INS ET RW 
Mm | Pl) — 2 5 Ei 70 een (5.52), 
where 
ß 2i d Fa 
26) = w&(o) = b 108042) +’ 7 | 
EDEN | a a ae 5.93) 
1 n ( 
ALLE 2 : ; | 
Des 
From equation (5.51) we have 
or Pin 
Am | vd) — res u ee (5.54). 
Equating expressions (9.52) and (5.54), we get 
ao Wie, | 
a YA 3 Eee ‚55). 
G=2[l0 5,20 (5.55) 
Substituting this value in equation (5.51), we have 
BEN 1/0 @\, © es 
Wu) = 3b) + Ars „)2@) ee (5.56). 
Substituting this equation in equation (5.41), we obtain 
en BDA NEOR Yrio; = © 
—) =, z) + — (— Z r LEW 
u An 2 2 4 2 Fr =) ©) 2 = 29 


For numerical calculation we shall express this equation in the form of trigonometric series 


re P | @ BERN ENTE TEN EN ns BI ncos2nrv ! (5.58). 


4n 20 Aw(l—cos2n)) Exp 2ngid)— 1 


The series in this expression converges when the absolute value of £, is less than unity, and 
as &, approaches zero it converges more rapidly. 

Particulary, the bending moment M, at a point u =2ma, and =2na, or at the 
center of a column is caleulated by putting 


IN) s j pP 
gt + log .) in place of 4-{F10,0) +51 Kı(0, 0)} 
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in equation (5.1). Similarly the bending moments M, and M, at the point are calculated by 
putting zero in place of ne {F;(0, 0) + &k K;(0, 0)}, where j=2 and 3 respectively.: Hence 
7 


we obtain | f 
P 2 90, 1 MR Ipbugägaun | 
[ 1 | 
= u |Rf2I0s Qu) 4°] + 2108 Ua} ka 
dr 6 2 
Mer Zu ZEIGE, SR). 
1 HM ze (5,60) 


6. Some Properties of the Functions B,(z) and Y,(z) 
In the preceding article we have derived the functions D,(z) and Po); or Dol&r, Es; © ©3] 
and Yl&ıs &; @1, @;]. Let us examine some of the properties of these functions. 
(i) From equations (5.27) and (5.51) we can easily show that 
exp{D(—- = —erp{DM ---: ren (6.1), 
exp {D,z +2 0,)} = — exp {Dulz) + ridz}, 
exp {Du(z + 2 0,)} = — exp {Dı@) —riEı} 
VID) WVlD" sea m re ee (6.3), 
er De (6.4), 
P& +20) = Pol) 
(ü) Let us denote by D,(&,&3; ©), ®;) and PE,&;; ©, ©;) the functions 
DE 835 0 ©) = Ri Do lEr, E35 © @s]} » 
— DE, &35 @1 ©) + I Dy(&, &35 @1 5) 


ji a‘ ß) ei n) u RE 
ee @ do, + @g a DulEı &; ©, ©3] 


1 
Pl; &; @1, @;) = et 
— Pole &; 1 @;) + I Pl, E35 Op @5) 
>= ol&ı, 3; @1 @;] ’% 


where Wl&, &3; ©], @,] is the function obtained by substituting 2&, ®, + 28, w, for zin equation 
(5.56). Then the bending moments can be expressed in the form: 


D > 
Mm= 7 DE, 85 0, OS) + RP a a (6.7), 


for ee 


is 


a7 | 
mn En 


i=e1b233; 
(ii) As the bending moments must be even functions of z or 2&, @, + 2&, w, we have 
D— &, — &5; 0, ®,) = Di 35 01, @y) , 
Pl— En — 3 0 9) = Plln &ii O1 Du) 
(iv) As the functions which represent the bending moments must have the periods 2 &, 
and 2 w,, we have 
DE + m, &, + n; 0,0) = DE, 3; @, @,) 
Fe +m,&+n; 0, 0) = PlEn &; 01.0) 
where m and n are positive or negative integers or zero. 
(v) As we see from equations (5.27), (5.37) and (5.56) that the funetions PD, and YP, consist 
of dimensionless terms, we have 
DE, 53;00,C@)= IE &3; 01 @;), 
L .p » . 
Pre 8; 0, co) = File 60: 0) 


where c is any real number. 
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(vi) In the preceding article we have calculated the bending moments by drawing the 
network of parallelograms with sides 2 w, and 2 w, on the slab. We have assumed that the slab 
is supported by columns located at the centers of the parallelograms. Now we shall cover the 
slab with another network of parallelograms with sides of magnitudes and directions 2 @,, and 


Br ee > 
Se 


Fig. 10 


2 @„, on condition that the center of each parallelogram is located on each column (see Fig. 10). 
As the centers of the latter parallelograms coincide with those of the former, the following equa- 
tions must be satisfied: 


o, = Mm, ©, + my, @;, (6.11) 
BR ee LH; 
0, = N, ®, +4 N, @y 
where m,, m,, n, and n, are positive or negative integers or zero, satisfying the equation 
m, m; | 
ER ee ge (6.12). 
A N, 


Any point z = 28, ®, + 28, ®, on the slab can be expressed in the form 


NE PEN RE N (6.13), 
where (27% 
m, m 
EN AD EI EEE ea (6.14). 
n, N; & E35 


As the bending moments of the flat slab are not influenced by the shapes of the parallelo- 
grams, provided that the columns supporting the slab remain in the same positions, we have 


Dyle1 &s5 ©1 0) = DrlEm En; Om On) i: NE (6.15). 
FE, &35 915 @;) = Diez n? Om» on) 
As a particular case, we have 
ne 6.0) 
PilE &35 01 @;) = Pl — 1: Or 0) 
(vi) We know the four types of the Theta-functions, and remembering the relations be- 


tween these functions, we have?) 
d,(0|r) 


1 
D,(z) = log 9,(v|r) 3 log > BE N (20 
1 (0 z i 
D,(z + o,) = log dy(v|r) 3 log ’ z = +ri&v-+ .n 
1 d(0 ; i/1 
De + 01 + 0) = log dw) — log nn +ai0+ Tr & + 5), (6.17). 
100) mi 


Doz + ©) = log dılulr) N log In +rid;v+ 2% A) 


Taking the real part of these equations and substituting zero for z, we have 


1 d(Olr 
D, (5 ’ 0; 915 oo) Ze RN j'og d,(0lr) 3 log n : = 
Be 1, 00) 
D, I; s 5 0 2 — N j'og d,(0|7) er 3 log Fr | ’ N (6.18). 
1 1 (Or 
2,(0,5 50104) = 810g AO) ze | 


; : 2 SER 2  ı; u : \ce ZAMM 2 
5) A. NAvaıtransformed his solution expressed in a trigonometric series into the d,-function (See ZAMN 
aa We can easily deduce NApar’s solution from the second of the set of equations (6.17), by substituting 


a purely imaginary number for 7 in it. 


21 
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/ 


Hence ug IR ’ 
0, (5,90 0) + 9 (9550 0) + 9(0:5 50% o)=10g2 . . . (6.19). 


(viii) Now we shall consider the relations between the bending moments at the points 
on the coordinate axes and those on the lines dividing the span of the flat slab into two equal parts. 

The bending moments at point A (see Fig. 11.) on the flat slab carrying a load of intensity 
PJ/S uniformly distributed over the entire surface of the plate, supported by columns indicated 
by the white circles in Fig. 11. and receiving reaction 
from each of them, are caleulated by using coordinate 
axes U,, u, and are represented by the expression 


e el E04 oo) +kY, (7 63; De) 


(6.20), 


M,; 


EHEN 


When the flat slab is supported by columns 
indicated by the black eircles in Fig. 11 instead of 
the white, the bending moments for the point A 
are calculated by using coordinate axes u, u, and 
are given by the expression 


| 
M; = 7 18;0,8,50, 0) +3 kY0,8501@9)} (6.21). 


Fig. 11 Finally, when the flat slab carrying a uniformly 

distributed load 2.P/S is supported by columns indi- 

cated by the white circles and those by the black circles and receives reaction P from each of them, 

the bending moments for the point A are calculated by using coordinate axes ız,, u, and given by 
the expression: 


pP [0) \ 
u 19 (0 6:5, 0) +20,k%, (0 &; 3 oo) ET (6.22). 


Asthe bending moments given by equation (6.22) must coincide with the sum of the bend- 
ing moments given by equations (6.20) and (6.21), we have 


1 1 
7 & E55 0 or) A KG ip: SER or) + ®,(0, &;; ©, @,) 
| . 16,23); 
=, X k 7,0, Es; 15, os) = ®; (9 SR - ’ or) + 2 X k y, (% Er - os) 


Hence 
[43] 


1 
Der = 1 = £ 
D, = » &95 On or) = B,; 0, Eos rw 0.) — D,(0, S35 LOFR ©) 


jr 


[42] \ 


1 
5 N) ._ 7,0, E35 01 ©3) 


| un 


Yo ‚S3; On o) =27; (9 $3; 


Particularly for &, = 0 


1 2 (012 7) 
9, (>, 0; 0, DE ren 
(go m) =Ig2+ "1° 30 |’ 
1 ; . 1 1 /o, 2) 
—®D, © BE PR oo) 7 95 05 | =. er ae r, 1202 DZ}, (6.25). 
— De E ‚0; 0, os) +1 e ‚0;o os) = —-. an -- I u a $ 12202 z 
2 b P ’ 13 a w © 2 ar 2) En (2 T) gt (r)} 


Using the similar method, we have 


1 f D: 
BD, [& I ‚ 15 oo) = D, (& D O1 u DilE,, U» 1, 3 ) 3 

| | | " us (6.26). 
P (& 300 or) =2W, (& 0; 0, 2° FE: 0; 0, @;) 


4 
f 
j 
}- 
y 
* 
4 
’ 
% 
j 
s 
. 
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Particularly for & = 0 


et 
2, (0, gl >) 
De LE TEE) = en 
( a ) a TE 
1 3 1 1999; — 
(0, 550% 0) +1 9(0,5 50, a As): (6.27). 
Du 17 u ee, T 
7, (0 DER oo) +i 7.0.5 01, 2 = >= = 325) -ze) 


(ix) Remembering the definitions of the bending moments M,, M, and M,, we can easily 
show that 
| DE 83; 01 e, @,e) = Bule,&35 01, @3) , | 
DslEr &3; 01 €, w; e'°) = Dy(E,&; 1, @) 608 25 — DylE,,&; @, w,)sin26, (6.28). 
DylE, &5; 01 €°, w; €) = DE ,&; @, @;)sin26 + DEE; , @,).cos 2 6 | 


Dr 2220 er) SUITE 2: 00,)008 20 — DE, 3.0, 0,5 20, 6.29) 
NEE ne 0 2,06) = Yl E00), 20 Pl: 0,0) 6052e 


where ö is any real number. 
(x) We can represent equations (6.11) and (6.14) in the form of modular transformations 


Walsh “ vi ee 
=(_ r N Ei Wr FE (6.31). 


Generally the modular transformation can be expressed as a product of the elementary 
transformations 


m m; [1 b OT 0... 1, /1.b, RN 
el En ar I 
where by, d4, - . . , b, are integers or zero. Then we have 
i 
exp [Del Em Eins Om On} = exp |BılEı 650, 0) — 8A + >) ss). 
where 
DIE ED er wen Ad A er en ERRRTER, (6.34). 


For the particular case where & = 0 and &, = 0, we can deduce the following equations 


i 
r en Es 
IQ, )P exp FE 2 —= (m, + my r) {[Q,r)P exp ir (— Yb) 2 10.35), 
or 
9,0) = (m, + m; 7) {H0|7)}? exp = (AEEN | BE . (6.36), 
where 
N SE (6.37) 


m + mt. 


Substituting expression (5.27) for the function ®, in equation (6.33) and by using for- 
mula (6.36), we have a modular transformation of the Theta-function 


AN ; 
9)? = Won)? (mı + mar) exp 2mzrivv + (— N | (23) 
where a 
_ — ee) 
30. er 


21* 
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Differentiating the logarithm of equation (6.35) with respect to r, we have 


ZI _ (m m mem tm) 2... 640) 


T 2 Ad 


(1 A 
For the elementary transformation _ 1 0) we have the wellknown relations 


1 (z = )' - ir{ö,(u)PP exp (2=i -) a aaa: &- (6.41), 
(0) = tion re 
Ze +2(—)+1= ER (6.43), 
-OrT 1 
a Sie.) 0, So re (6.44). 


2 


7. Bending Moments in a Flat Slab with Rectangular Panels 


We shall now discuss the bending of a flat slab with recetangular panels. In the present 
case the oblique coordinates u, and u, become the rectangular coordinates, and we take the r- 
and the y-axes as they coincide with the u, and the u, axes respectively. Then we have 


yı=d; re 


oo, =dq4;, VO, —=1d,, 


t=U, y=U, 


As mentioned in the preceding articles the bending moments can be expressed by the use of 


& &, @, @,; and k. Particularly for y,=0 and y, -5 we shall use notations &, n, fand k, 
where 
E = & $) N =< & ’ 
De ee EEE . . . . . . . . . . . . (7.2). 
3 4a,d, 


Then the bending moments in the flat slab are expressed as the functions of £, n,tandk. 
We must note here that they are also functions of a, and a,, but we shall omit these letters assum- 
ing that they have always constant values. Then we can represent them as 


P i er . 
M, = a U(&, nt k) [= 128 3 I Ss 
If & and 7 satisfy the condition 
IE— MD +LÄNn— MD... 2... 2 20. (7.4), 


for all the values of m and n, we can divide M, into two parts, one of which is independent of 
% k and the other has the coefficient a, k 


Pp N j 
M,= Te {WE n5,0)+0,k Dam d} won aan va 
Comparing equations (7.5) with equation (5.1), we have 


0459560) = Bd; 1,it),) 
Ve) re Lid |  e 


As mentioned in the preceding article, we have 


Ul&E +m,n+n;t, k) = UlE,n;t,k), Tel 2 3 Deere (7.7). 
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ag equation (7.1) in equations (5.2) and (5.3), we see that U;(&,n;t, k) has the 
following properties: 


UY-&1t = Uen;th), ve 
ee ne = 8 
U—&,n;t, k) = — A k), 
neuen EN N +. (7.9), 
1 \ 
U; N; % k) = U (m 38 Fr h), 
1 
Une, (m :-, k), een (7.10). 


U& m. = Urn) 


For the points which satisfy the condition expressed by equation (7.4), we can calculate 
the bending moments by using the following expressions: 


1 \ 
Em) =-—atlg+ m) 
+5 1og {2 (eosh 2a 17 — cos 278)} 
De ru 
32 n{exp(2nat)—1} ° 
1 atnsinh2rntn 
5 u wi 2 
6,9; 40) atlgt m) cosh 2% in— cos2nE 
cos2nz&sinh2natn 
Aa 7-11), 
Hain) exp (nat) —1 > 
ehr 
ey, 
aıy cosh2nnt—1 ; 
’ » atnsin2n& 
Ur, 936,0) cosh2rtn—cos2n& 
=, sin2nn&cosh2nztn 
Be "2 Senna 1 
sin2nnz&sinh2nntn 
er a1) cosh2nat—1 
e 1 
VEnd=- 7: 
Ä 1 cos2r&cosh2rtn—1 
BEI SEFTR och min cos InER 
ncos2nz&cosh2nntn 
Hal) en Te (7.12) 


sin 2rz&sinh2rtn 


VE —n (cosh an In— cos 2 &)? 


n sin2nzd sinh2natn 
Ham 3" exp Ona)—1 


The series in these expressions converge, when the absolute value of n is less than unity. 


Particularly for &= 0 and n = 0 we have 


rs a u ee 
£ IE nt 
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2 l il —. k— 0, 
0,0, 054,1 = log4akt + 210g Ge) G—Z 
1 £ Art: 
0,0, 05, = —(l—ok) 2; + Zi 9) 


U.0,0:5, =) 
From equations (6.24) and (6.26) we get 


Ujlg L, R) — U;(0, 1; 21,2) — U0,n;t, k), 


U,(6.5; 1, k)= Uj(& 0; = 2) — Ud&,0;1,k), 


1=1,237: so 


Using the similar method and remembering the properties expressed by equations m 
and (7.10), we can deduce several relations which exist between the bending moments produce 
along sections dividing the span of the flat slab into three, four and six equal parts. We must 
note here that the following relations are applicable only for U, and U,, but not for U,. 


il 
ug m6R)=5 Ulg-n 2628), 
’ 1 Nr eh Be (7.15), 
(68456), Ul&555>2%) 
1 
Ulzn; L, =, U;(0,531,3K) —— Uj(0,n;t,k), 
1 1 ! er A 
Ul&gs6hl=5 U(505 3,34), VE 06h), 
ae, a je 1,2% Se 
U,(gn:6k) 5 Ul5, 73438) Uılg meh), 
a 1 \ 1 et 1 j 
ul sk), Ulbg5 glg Ul& sth), 


Particularly for {= 1 we can deduce the following relations: 


ee N Seen Dane 
un, = Ui > E) tl 9 en 2 


Ee—n E+n k te +- 7 Im TE 
Uz(&,n; L, Ki ul 5) 3 9 sl; SE: ul p) > 9 BR a: 
EN ee UV 1679 1.& 
a ea ee va 
Using the above-mentioned equations, we see that the bending moments M,, M, and M, 
en 1. 
64° 90.2 u No 
can be deduced from the bending moments at the point &= 0 and n = 0 of flat slabs of which 


ee 
Br 4° 3: 5 ’ 1, 2, 3, 4, 6. 


(7.17). 


- 


1—E—ı k 
Use, N; T. k) Hui U — Ba . ) 


of a flat slab with square panels or {= 1 at points& —=(, 


8. Numerical Calculation 


As we encounter most often flat slabs with square panels namely {= 1, the values of 
U;(@&,n;1,0) and U;(&,n;t) fort=1 are given by expressions (8.1), (8.2) and by Tables 3—7. 


U,(0, 0; 1, k) = 0.738 168 + - log k — 7 k—o,, 


Us(0,0;1,)=0, ERBE 
U,(0, 0;1,K) = 0 

‚ 1 

Un; )=—— 7, es u BO 


2 
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Table 3. Numerical values of U,(&,n;1,0) 


— ._ u en, 


—_—_—_—_—_ä_ä_—_—__——— 

£ a En 2 ER 1 5 DR 
a E=0 | = | =; | = | i=z | = =; 
a ne a 7 


—1.185320 | —0.525468 | —0.177022 0.027580 0.138119 0.173287 
—1.185320 | —0.849465 | —0.423930 —0.133218 | 0.050715 0.152987 0.185878 


—0.525468 | —0.423930 | —0.219493 | —0.029186 0.109746 0.192249 0.219493 
—0.177022 | —0.133218 | —0.029186 0.086643 | 0.182173 0.242964 0.263665 
0.027580 0.050715 0.109746 0.182173 | 0.247073 0.290712 0.305975 


| 
0.138119 | 0.152987 0.192249 0.242964 | 0.290712 0.323997 0.335855 


0.173287 0.185878 0.219493 0.263665 | 0.305975 0.335855 0.346574 


Table 4. Numerical Values of U,(&,n; 1,0) 


0 | =. | E= 4 = | 8=2 | E=2 | E=5 
0.503178 0.511843 0.523643 0.535398 | 0.543976 0.547110 
—0.503178 0 0.308534 0.419693 0.470897 0.496223 0.504120 
—0.511843 | —0.308 534 0 0.202235 0.315874 0.373953 0.391939 
—0.523643 | —0.419693 | —0.202235 0 0.142283 0.223724 0.250088 
—0.535398 | —0.470897 | —0.315874 | —0.142283 0 0.089604 | 0.119905 
—0.543976 : —0.496223 | —0.373953 | —0.223 724 | —0.089 604 0 0.031190 
—0.547110 | —0,504120 | —0.391939 | —0.250088 | —0.119905 | —0.031190 0 


Table 5. Numerical Values of U;(&,n; 1,0) 


ma | era | et | et | et | 5 | 
\ — 

0 | 0 0 0 0 0 
0) 0.471987 0.344 842 0.218930 0.129028 0.059952 0 
0 | 0.344842 | 0.391939 0.304033 0.195969 0.094696 0 
0 | 0.218930 | 0.304033 0.273555 0.192189 0.096 944 0 
0 0.129028 0.195969 0.192189 0.143459 0.074851 0 
0 0.059952 0.094696 | 0.096944 0.074851 0.039856 0 
0 | 0 0 0 0) 0 0 
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Table 6. Numerical Values of U5(&,n; 1) 


Manuskripteingang: 28. 9. 1960 


ESIESIEC ESSEN 

N = 0 — 22.928762 | —5.771474 road _ 1094220 

n= 5 22.928762 0 — 2.781357 | —1.916483 | —1.346921 —1.074374 | —0.994539 

A . 5.771474 2.781357 0 —0.727106 , —0.801024 | —0.764668 145437 

Me — 2.641680 1.916483 0.727106 0 —0.312197 | —0.425609 | —0.453 241 

m = 1.606194 1.346921 0.801024 0.312197 0 —0.159513 | —0.207455 

N n 1.203892 1.074374 0.764.668 0.425609 | 0.159513 0 —0.052219 

n= = 1.094220 0.994539 0.745437 | 0.453241 0.207455 0.052219 

Table 7. Numerical Values of U;(&,n; 1) 
en: eG | = | E=2 | e=4 | e=+ |, E25 

ven 0 0 0 | {) 0 

n= = 0 —11.438272 | —3.625066 | — 1.310922 | —0.542967 | —0.207455 

Ye = 0 —3.625066 | —2.782010 —1.503283 | —0.745437 | —0.308252 

" = = 0 —1.310922 | —1.503283 _ —1.094220 | —0.645465 | —0.291127 | 0 

1, + | 0 —0.542967 | —0.745437 | —0.645 465 | —0.430 378 | —0.207 901 | 0) 

= > 0 0.207455 —0.308252 | —0.291127 | —0.207 901 | —0.104676 | 0 | 

4 | 

72 : g 0 0 a u | 
j 
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Von Kärmänstraßen induzierter Flüssigkeitsttanspott 


dei Hir* Von Br. DOLAPTSCHIEW u. Iw. TSCHOBANOW 


E EN Mit Hilfe der Formel von Lichtenstein wird der von einer Kärmänstraße induzierte Flüssigkeits- - 
- transport durch eine gegen ihre Achse beliebig geneigte und zur Stromebene senkrechte Wandstrecke bestimmt, en: 
in bezug auf ein Koordinatensystem, welches mit der als ruhend vorausgesetzten Flüssigkeit verbunden ist. Das ” 
Problem wird auf eine Integralfunktion zurückgeführt, welche eine Verallgemeinerung der Funktion von 
Lobatschewsky dargestellt und deren Reihenentwicklung nur für den Fall gilt, daß irgendein Ende der 
Wandstrecke auf einer der beiden Wirbelreihen liegt. Für die übrigen Fälle wird die Fouriersche Ent- 
h wicklung und die Entwicklung der .allgemeineren Integralfunktion in eine Potenzreihe gefunden. Es werden 
4 effektive Formeln für den gesuchten Flüssigkeitstransport angegeben. 


m Lichtenstein’s formula is used to determine the transport of liquid, induced by a Kärman street, through, 
i a wall normal to the plane of flow and tilted by an arbitrary angle with respect to the axis, the system of 
DE reference used being the rest system of the liquid. The problem is reduced to an integral function which is Be 
ei a generalized Lobatschewsky function with a series representation valid only in the case where one of- ü 
; the ends of the wall coincides with one of the two vortex lines. For the remaining cases the Fourier series 
E. and the representation by a power series is found. Effective formulae are given for the transport of liquid. 


ABTOPbI HAXONAT TPAHCHOPT ;RUNKOCTH, UHNYIMPOBAHHbIuH KapMaHckofi BUxpeRoi NOPosKk- 
Ba Koi yepe3 BOO6PaskaemyIo CTEHY, 3ARJIIOYAIOINyIO IIPOM3BOABHLEIH YTOJ C OCbIO MNOPOFRKU 
| MH HEPIEeHAHKYAAPHYIO K INIOCKOCTNU TeyeHua. us 9Toi Wenn Mcenonbsyerca PopMmyaa 
JInxTenumrteiHa OTHOCHTeNBHO CPABHHTEeIbBHOH CHCTeMbI, CBA3AHHOÜU C >RHIKOCTbE. 
3apaya lIpuHBOAHT K HHTETPAANAaM, NOAMHTErPAABHAA PYHRIMA KOTOPbIX ABJIAeTCH 3HAyH- 
TeJIbBHbEIM O0600IMeHHeM PYyHRUMM Jlo6ayueBcKoro, PA3BbITMe KOTOPOoÜU OTHOCHTCA TOABRO 
R CJAIy4yal0O, KOTAaA ONAHH KOHEIN OTPe3KAa-CTeHLI HAXOAHUTCA HA OMHOH H3 BUXPEBEIX IeNOyerR. 
Ist OCTAJIbBHEIX CIIyyaeB aHo PasBuTnme B pn Dypme uU B cTeNeHHBIi PA 060Ö60ImMeHHoH 
#4 pyarımm. Haüpensı IhpPekTuBHBIe PoPMYABI HIA HCKOMOTO TPAHCHOPTA FRUAKOCTN. 


Eine zweiseitig unendliche Wirbelkonfiguration, d. h. zwei parallele Reihen von Wirbeln in 
4 symmetrischer, in schachbrettartiger oder in asymmetrischer Anordnung, mit der Straßenbreite 
2 h, der Wirbelteilung 2 ! und der Reihenverschiebung 2 d, mit der konstanten Zirkulation — /', 
bzw. + /’und der konstanten Translationsgeschwindigkeit U oder W(U, V) der gesamten Wir- 
belkonfiguration nennt man — wie bekannt — KArmänstraße — ein Modell der physikalischen 
Erscheinung eines schlecht umströmten Körpers (Bild 1). 


ED IN, 5 2! Ken 
2h D— Eh Ü— 2h NS 


ER symmetrische schachbrettartigeAnordnung er 
Anordnung 
Bildı = 


Dieses Modell ist bis heute zweierlei theoretischen Untersuchungen unterworfen worden: 
Einerseits wurde das Problem der Stabilität der KArmänstraße vielmals diskutiert, anderseits 
das Problem des Widerstandes, den der umflossene Körper erfährt, ausführlich behandelt. 


i Man kennt die von KArmAr [1] gegebene Stabilitätsbedingung 
l 
sh(#m)=1, = 7» 0,281 a IE 


für eine schachbrettartig angeordnete Wirbelstraße oder die von MAve [2] und Dor.aptscHıew [3, 4] 
gegebene allgemeinere Bedingung 


sh(&r) =sin (An), U . (2) 


d 
I 
für die asymmetrische (zweiparametrige) Wirbelkonfiguration, die für A = 0,5 (d = 1/2) in die 
KArmäAnsche Bedingung übergeht. Auf diese Bedingung gründet sich auch die KArmäÄnsche Wi- 
derstandsformel. 

Später sind aber neuere Untersuchungen von KAMENKoWw [5], KorscHin [6], CODDInGToN [7] 
u. a. erschienen, die sowohl die Stabilitätszahl als auch die Stabilität selbst und daher auch die 
Widerstandsbestimmung angezweifelt haben. 
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So hat z.B. Kamenkow 1935 die Wirbelstraße kleinsten Widerstandes mit der Zahl 
h/l = 0,245 statt mit der KArmAnzahl h/l = 0,281 berechnet. Ko 

Andererseits zeigte KorscHin [6] auf Grund der Lsarusowschen Stabilitätstheorie, daß 
bei Störungen höherer als vierter Ordnung überhaupt keine Stabilität existiert (DOLAPTSCHIEW [8]; 
[9] hat früher bewiesen, daß bei Störungen zweiter Ordnung eine solche Stabilität immer noch 
vorhanden ist) und daß die Stabilität KArmAns bei Störungen erster Ordnung nur die kleinste 
Instabilität bedeutet. . 

So sprach das Modell einer zweiseitig unendlichen aus unendlich dünnen Wirbeln konstru- 
ierten Wirbelkonfiguration. 

Dieses Modell mußte allerdings noch von anderen Seiten her geprüft werden, sowohl theo- 
retisch als auch experimentell. en 

Am IX. Internationalen Kongreß für angewandte Mechanik in Brüssel haben wir über einen 
anderen Effekt dieses Modells vorgetragen [10], [11], nämlich über die von der Wirbelstraße auf 
die umgebende Flüssigkeit induzierte Bewegung, welche dem Seiche-Effekt CHryYsTL-HALms [12] 
ähnelt, bei dem das Niveau des Genfersees periodisch pendelt. 

Unser vorliegender Aufsatz betrifft eine andere theoretische Frage, die ebenfalls mit dem 
betrachteten Modell einer Wirbelstraße im Zusammenhang steht, nämlich um die Berechnung des 
Flüssigkeitstransports, welcher von 
einer KArmänstraße durch eine be- 
liebige Fläche S, während der Zeit 
{, — t, induziert wird. Die Auffin- 
dung dieses Flüssigkeitstransportes 
verlangt die Verwendung der Defini- 
tionsformel (nach LICHTENSTEIN 


13) | 
Q = [at ff et, 7) w,(t, 7) do (8). 


Hier bedeutet o die Dichte und w, 
die zur Wandfläche S normale Ge- 
schwindigkeitskomponente desjeni- 
Bild 2 gen Flüssigkeitsteilchens r, welches 

sich zur Zeit? auf S befindet (Bild 2). 

Wir beziehen die Betrachtungen auf ein Koordinatensystem o x y = [o], welches mit der 

in Ruhe vorausgesetzten Flüssigkeit verbunden ist, und betrachten zunächst die ebene Wand- 


ee = A,(&) As(&;) mit der Höhe eins, die mit der Wirbelstraßenachse den Winkel x = 5 
schließt. 


BL Die Flüssigkeitsmenge der als homogen und inkompressibel vorausgesetzten Flüssigkeit 
ist dann 


t 


dy . de 
02, 2)= j dt [eos sn 7) dr: on VS 
wobei 5 
de RE dy d 
IE == he 7 lz) dt = Im 7 f(@) Pe 2 (5) 


die Geschwindigkeitskomponenten eines Flüssigkeitsteilchens sind, welches nicht mit irgend- 
einem Wirbelzentrum zusammenfällt, und 


RE; I ! 
pi, 95 [? le aan, h+-U ) 
5 e 
> 


i 
sin z+— +ih+U 
tg Hit) 
das komplexe Strömungspotential in bezug auf das System [o] ist. 
Man erhält hieraus für die Geschwindigkeitskomponenten die Ausdrücke: 
IT 
a na ly+ Den 
LERNTEN A en 
d Al "| 7 FR 
kr ch — [y+ — D*h]— cos — (e + Ut 
I I ( (7). 


dk TV 
— =-— sin - 


TAT TER DIT > 
it! h7y+C En—coT@+ Ui) 
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Mit den Ansätzen 
SW Ta+UD, dr 
0.) = CI) “ (=. BR), re er. (8) 
= chT [tgaası +m+-1Eh], (Re) 


nimmt der Ausdruck (4) für den Flüssigkeitstransport Q die Gestalt an: 


5 9,52) 
Fol a 
> I = > () | 1n.(d; 2 Sin 0;,) de. (9). 
Kl gl 
9;,() 


Zur Berechnung der Integrale in Formel (9) ist es im Falle a;,; > 1 zweckmäßig, die Funktion 
F(&) = [In (a + sin a) de, Pr A (10) 
0 


zu betrachten. Diese Funktion stellt eine Verallgemeinerung der bekannten Funktion von Lo- 
BATSCHEWSKY dar, nämlich 


” 


rn k 

Li) = — | moszde-am24 5 Se 1) 
0 

die nur im Falle a = 1 verwendbar ist, d. h. wenn ein Ende der Wandstrecke, z. B. A,, auf einer 


der Wirbelreihen liegt (Bild 2). 
Die Integralfunktion in (10) ist in eine FourIzzsche Reihe entwickelbar mit den Koffizienten 


Sn 


em: 


2n 


1 costsinn Tı 1 al! 
= de = — ger. 1 al ; 
”«, a8 Fran GE au Ewa) 
% PR 
1 [ coszcosnx 7 nr 1, 
— 2 Ze ee 2 Se 1 en 
zb, = ac): eiMııcıy 
0 
[ xcos« a+y®—1 
— = Er o). 
Ro=2nlina er —=2rlin - 5 (13) 


Die Auswertung der Koeffizienten a„, b„, (n = 1,2,...) erfolgt folgendermaßen: 
Nach der Definition hat man 


27 2n 
Duke Fer d. L fecos zsın NE 1 ] 
a n(a-+sinr) cosntdı = eure & 
0 0 
a =E el), 
b„= | In(a+sinz)sinnzde = SER dx 
Der z a Eee 
0 0 
m—l1,2, ) 
Es gilt aber [14] 
ee Be = era 
(rer Be Lane] oa, 
| costsinn® ED, Vaa— 1 (15); 
+ sin x ig) (n=2k+1) 
27 
" eosTcosnt 5, 0 (n=2k), 
i ZB — le 16). 
| a+sinz Ey (Ya — 1 = E Y@—1)°] eu (16) 
$ z 
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Mit | | 
a=ch$;, = infa +17 y@—1], K=12) ....: (17) 


gelangt man sofort auf die Formel (12). 
Die Auswertung des Koeffizienten a, ist wesentlich komplizierter. Aus (13) folgt 


27 
#0 = 2nna— | 04 an de el ee (18). 
Ö 


a+ sinz 


Zur Berechnung des Integrals in (18) betrachten wir die Funktion 


x zlogz - 562 
ur aa pe ei 


die bei a> 1 holomorph im Einheitskreise ist, aus- 
schließlich des einfachen Poles z, und der Strecke 
—1<sz=0. > 

Die Funktion (19) wird längs der Kontur von 
Bild 3 integriert; hier sind die Kreise |z| = 1 und 
|z| = e gezeichnet, wobei 0 <e< ||. 

Indem man setzt 


T=[Gd)dz, L,=/fÖle)dz, (i=1,2,3,4) (%), 


bei 


0] — A,As ’ 03 — A,ABCA, Pi 
sn O3 = A,A, ’ 04 — A,FEDA, ’ 
erhält man, nach gewissen Modulabschätzungen und anderen Maßnahmen, daß 


1 
i 2 x de 
ee Ar = —4inK 
ae » ae net 
0 * 
ET ee e1). 
lim I, = = 
I 2 | a-+ sin Ener 
lim I, — 0 j 
e>0 # 
Wir haben folglich ü 
; i .2, (In || +iargz,) 
EA K-UgiR BERN, 1 
a nn di 22, +2ia te a ae (22). : 
H 


Es ergibt sich 


1 1 
1 d(x? — 1)? x de 
K= al y Na 2 re 9; i : 
2 | @— IP +4IR(a@ 1) +4 F ia | (@— 1% FI =M +2iaN (23), 
0 0 


wobei 


re a — 1 
Me 5; Ir In (a — Y@—1) + 5 In2a, 
EIER 24). 
zb 1 N — — 
N= 7 Var = (7 y@—1-—-aarctg Va? — N) 
Aus (22), (23) und (24) berechnet man den Wert von L und daraus die Integrale 
"OsinOdd 2ma im ee | 
See ei E — arctg Ya — N) ne A (25), 
; 0 cos 0 dO ——— 
| at sind =2rin2a(a— Ya —1) ne EEE u - (26). 


—n 


Aus dem letzteren Integralwert ergibt sich die Formel (13) für a,. 
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Also ist die gesuchte FourIzgsche Reihe die folgende 


| u 
ar ur Fun a 


In (a + sinx) = In +2 


- r=]1 
| 2 A — 
x | +E-1WIeD° wenc+ll+-Wrtll-1) ? sinn ) en) 
Sie konvergiert für jedes x gleichmäßig und führt zur Entwicklung 
x 
@—1 
ee Re ! 
n-+1 
Nr De e "er| In Ce ale- rt en, Et Ir ie st) a cosn x) (28). 
n= EL 
Hier stellt die Konstante 
ee 16 
2 G(9.) = 22, @n+i ee (29) 
eine Verallgemeinerung der Konstante von KATALAN 
Nah: 
G N ee ten 30 
SS, @n+1)% ( ) 
dar. 


Anmerkung. Es sei hier bemerkt, daß die Funktion F(x) auch die folgende Entwicklung 
in eine Potenzreihe gestattet: 


% 


[ma + sin daran + 3,4 De KR (31), 
wobei v 
ra d’(u v) 35 
2 ee ee ae (32), 
mit 
dr(u v (7 IH a o 6 : 
an ei — 1)r > Se er—20) 
(33), 
und mit 
der if D) 
der+i = 
lan —q ki a 
un) er Dei a )ar—20+ perH: 
% (34). 


Die Berechnungen sind auch hier ziemlich kompliziert. 


A; 


nal fn 
ng za Be 
LEN = TI ee 


a) b) 


Bild 4 


Es gibt_verschiedene Möglichkeiten für die Lagen der Wandstrecken A,A,' bezüglich\der 
Wirbelkonfiguration, welche auf die folgenden wesentlich | verschiedenen Fälle zurückgeführt 


werden können (Bild 4): 
y„>h, h) —h<y<h, cd) Ge, d —h<y<h, 


35). 
%>h; „>h; „>h; —h<y<h a 


j 
{ 

I’ 
h’ 
S 
Bu 
er 
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Diese Fälle haben folgende Bedeutung: es läuft in manchen Fällen durch die Strecke A,A, in 
einem Moment ein Wirbel, d. h. dann verschwindet das Argument des Logarithmus in (9), und 
das Integral wächst unbeschränkt. Die hydrodynamische Deutung dieser Tatsache bei dem an- 
genommenen Modell einer KArmäAnstraße ist die, daß in der Nähe eines Wirbelzentrums die 
Geschwindigkeiten der Flüssigkeitsteilchen unbeschränkt anwachsen. Die physikalische Realität 
fällt aber nicht mit dem ausgewählten Modell zusammen: infolgedessen nimmt man für den 
Flüssigkeitstransport definitionsgemäß den Hauptwert des betrachteten Integrals an, welches 
im Sinne des Hauptwertes konvergiert. Liegt aber, wie erwähnt, irgendein Ende A, oder A, auf 
einer Wirbelreihe, so ist dann a;; = 1 und die Integrale in der Formel (9) sind nicht von der 
Form (10), sondern von der Form (11), in welchem Fall die Funktion L(x) von LOBATSCHEWSKY 
zu benutzen ist. 2 

Man hat endgültig aus (9) und (28) für den gesuchten Flüssigkeitstransport die effektive 
Formel erhalten: 


To x ei ent Dipl ir Ba 
RE) are +kjn (a, + Yar—1) 7m 


n=1 


2 2 n 
x _1Y+r 6 "| If1 0 we Ba U 
Zee item Do a+z6+h 
1 1 n+1 Zi k cı nz | U 36 
++ ED yE eYanlutzatnl + @0). 
Für den Fall a) in Bild 4 ergibt sich die Flüssigkeitsmenge 


Doll Sn bes ana U 
LRZE %) = mu n? sın 91 (k —h) 


2 —nay; n - 
k 2, ar ji er 008 7 + E (h+ A) a 


i=1 


n—1 


TR IE 
+1+e- Pt) 7 sin B - Su + A) ch = IF 
Für den Fall b) hat man 


0 TolsS1l na U 
LOSE x) = 91 (t, E t,) (h —— y,) 4 U } n? sin 21 (la t,) 
n= 


—nnh 
nz, j Nzt 


x I +1] 1)? I- ee am. B + = (, + A) 


—Nnay 


ua nsch nat 


J1 U 
+ e sh gr cos = I + 5 (, + A 


n—1l —_nnah f 
ak. Ar u un Ei nzy . ma U 
1 + — DD" r!](—1) | u Be 1 Iı sin T I aan (+ on) 
R u ı nz h . nz | U 
ch — in %+5 (+ in) ei - 


en ausgeschlossene l“all, bei dem x = /2 ist [14], führt ebenso zu der Reihe (28). 
oe ne des Transports durch eine krumme zylindrische Fläche S mit den Enden 
A a ührt zu der Berechnung der Größe Q nach Formel (36), sobald die Bogenstrecke A,A 
eine solche Struktur hat, daß das Integral (6) einen Sinn hat, [15]. 47 
Ei Se ae h ransports bezüglich eines Koordinatensystems OXY = [O0], das mit 

m Beobachter, d. h. mit dem umströmten Hindernis ver N i \ 

| | s ernis verbunden ist, wird auf d 

Weg ohne weiteres geleistet. Daun 
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Fehlerschranken bei Eigenwertproblemen 
normalisierbarer Matrizen 


Von JoACHIM A. NITSCHE \ 


Die Vektoren x; mögen die Bigenvektoren einer normalisierbaren Matrix A approximieren und 4i 
die Bigenwerte. Mit Hilfe der Abweichungen |A ti — Ati und des Mindestabstandes der }; werden Fehler- 
schranken für die Genauigkeit von x; und A; hergeleitet, zu denen für x; nach passender Normierung noch 
ein Faktor Y3 FI hinzutritt, der durch die Nichtorthogonalität der Eigenvektoren bedingt ist. 


Let x; and A; be a set of approximate eigenvectors and eigenvalues of a normalizable matrix A. With the 
aid of JA — Arzı| and the minimum distance of the A; bounds for the error of x; and 4, are derived, those 
for xi being supplemented, after suitable normalisation, by the factor Y3 \ri? resulting from the non-ortho- 
gonality of the eigenvectors. 


IlycTb BeKTopb % U umcaa A; CyTb IPHÖJIMKEHNUA COOTBETCTBEHHO COÖCTBEHHEIX BER- 
TOPOB U COÖCTBEHHLIX 3HAYeHHI HOPMaJlusyeMoä MaTpuusl. C HOMOMBIO OTKAOHEeHNÄ |A Li — 
— A; | M MUHHMANILHOTO Pacroamnn uucest 7; BbIBONATCH TPAaHHIBI IOTpelIHocTeh IH TOU- 
HOCTU T; U A, HPMY&M LA %; MOocAHe MONXONAMeli HOPMUPOBKU N00aBJIAeTCH MHOFKUTEIL 


VE 2i|?, koroppiii 06ycHaBJImBaeTchH HEeOPTOTOHAALHOCTLIO COÖCTBEHHbIX BEKTOPOB. 


In einer früheren Arbeit (ZAMM 39 (1959), S. 322—325) wurden Fehlerschranken für appro- 


ximativ bestimmte Lösungen des Eigenwertproblems einer symmetrischen Matrix hergeleitet. In 
der vorliegenden Note sollen entsprechende Abschätzungen für den Fall nicht-symmetrischer 
normalisierbarer Matrizen gebracht werden.!) Dabei seien als wesentliche Punkte hervorgehoben, 
daß sich die Fehlerabschätzungen auch auf die Eigenvektoren beziehen und die Bauart der 
Schranken möglichst einfach ist. — Die gewonnenen Ergebnisse lassen sich in folgender Weise 
beschreiben: 


Die n linear unabhängigen Vektoren {y;} seien Approximationen der Eigenvektoren {X;} 
der n-reihigen quadratischen Matrix A. Die Zahlen {A;} seien Schätzwerte der Eigenwerte 
1A. — Die Normierung sei so, daß die Vektoren y; des inversen, durch 


yh=da (k-1L220..m. So 
definierten Vektor-n-Tupels Einheitsvektoren sind?): 
ein -i (= 12%,...30 Wer 
Weiterhin seien e; die Fehler der Eigenwertgleichungen von r;: 


s=lAn—Aul 22 


und es werde gesetzt 


s=Maug Wi 
g.= Min |, —Au| „0 u Wear 
ik 
Gilt dann die Voraussetzung 
4 
4neMax[{l, : ee 
VE x IE | ee A 


so sind die Eigenwerte von A sicher alle einfach, und es bestehen die Abschätzungen 


9 
A;—A;| £2neMax (7) TER ETW 
- (Al). 


Bei passender Normierung der Eigenvektoren {X;} bestehen weiterhin die Ungleichungen 


2 
KR: — il = 2/n’.eMax(1, -)- N@) a ne 
mit ; 


No=yzlıl. ni 


1 Re 7 LIpR » . n .. . ”. . ? 
) Es ist beabsichtigt, die folgenden Ausführungen in Kürze durch eine Kleine Mitteilung mit Zahlen- 


beispielen zu ergänzen. 


2 „rk N EN. Bla j" ® . 
) %* bedeutet Transposition und Übergang zu den konjugiert komplexen Komponenten von x 


re rn (DE 


Rn= 
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Es ist klar, daß die Fehlerschranken im vorliegenden Fall schlechter ausfallen als bei sym- 
metrischen Matrizen. Das kommt in den veränderten Potenzen der Dimensionszahl n zum Aus- 
druck und in dem Glied N(x), welches wegen der nicht gegebenen Orthogonalität der Vektoren 
{g} natürlich erscheint. Ergänzend sei angemerkt, daß die Zahlenfaktoren nicht die bestmög- 
lichen sind. Auch könnte grundsätzlich in den obigen Formeln n durch n— 1 ersetzt werden. 
Aus Gründen der Einfachheit in der Herleitung wie auch in den Endformeln wurde jedoch auf 
solche Verbesserungen verzichtet. 


ih 
Wir gehen von den zweimal n Vektoren {r;} und {y;} aus, für welche die Beziehungen 
Au-khtteisi,. st BE OR oe 20 (7) 
sowie (1) und (2) bestehen mögen. 
Durch den Ansatz?) En 
Sen dd LS, 


soll ein Eigenvektor von X gefunden werden. Ist A der zugehörige Eigenwert, so muß die Glei- 
chung E : 
AX—AX=0 


bestehen, die mit der Darstellung (8) in 
Au 2 ee A) Kay 3 ar &rök ii 


übergeht. Multiplikation mit 4; für = i und! + iliefert wegen der Orthogonalitätsbedingungen 
(1) für die n Unbekannten 
e=41—4; Eee re 


und &,...%—13 & 41» +. +, &, dien Bedingungsgleichungen 
0=&Mi tr ZPO ee 
K—hAr)n=EMi + Mir Een 
Hierbei ist zur Abkürzung 


EN Tee ER) 
gesetzt. Für später merken wir an, daß wegen der Normierung von }) und 3 (vgl. (2) und (7,)) 
Na BRUNO Br NE FRRRSR 12) 


gilt. 

= Sind Lösungen o und «, des Systems (10), das wegen des Auftretens von oin den Gleichungen 
der zweiten Zeile nichtlinear ist, bekannt, so stellt der durch (8) definierte Vektor X einen Eigen- 
vektor von A zum Eigenwert 1=4, + o dar. au Red 


> 


ie 


Wir wollen uns der Frage zuwenden, unter welchen Voraussetzungen über die- Fehler- 
größen {e;} in einer geeigneten Umgebung des Vektors r; ein Eigenvektor X liegt. Dazu müssen 
wir untersuchen, wann das obige System (10) Lösungen o, a; besitzt, die klein sind, wobei ihre 
Größe natürlich abzuschätzen ist. 

Unter Heranziehung der in der Einleitung definierten Größen e und q beweisen wir nun den 
folgenden Hilfssatz: 

Hilfssatz 1: Gilt die Voraussetzung (V), so besitzt das System (10) mindestens eine Lösung, 
für welche die Ungleichungen 


ol; Max |;| &2n e Max (1, ar 19 PB ae 6), 
Ki q, 
bestehen. ak, . 
Zum Beweise bedienen wir uns suczessiver Approximationen, indem wir eine Lösung von 
(10) gemäß der Iterationsvorschrift 


oa =-0 
tt = eu; + DI Wirerod We a RU) 
ke yarl= gun + I un enoh 


3) I’ soll die Summation über k mit Ausnahme %k = ı bedeuten. 


td 
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» 
zu bestimmen suchen. Wir wollen zeigen, daß die Summe }; c* e’ mit 
1 


für 0° und o? eine Majorante ist, d. h. genauer, daß die Abschätzungen 
ler — e!|, —or!lsee G(=12: 1, ee 
gelten. — Wegen |o!| < e || < e und c > 1 folgt die Ungleichung für 0! bzw. 0! — g° unmittel- 
bar. Weiter ist Bi 
1 


1; —4} 7 ot SU>TE 27 
wobei der letzte Teil aus 
< Ber: Penn. < ee. < ee 
€ =7Max(1,4j@) 4Max(1,27) 8 
folgt. Damit wird aber 
I 
of] S—:E Sce. 
„2 
Wir wenden uns den Differenzen 
vr — oP _— Pa 
WR Pr a re sah 


mn 


zweier konsekutiver Iterationsschritte zu. Diese genügen den Gleichungen 
BEER — N up w? 
K—A+etr)wtri= 3 ur wg —of vr+ 
woraus die Abschätzungen folgen: 


or+!l<n.ewer 
(q — joe +!) we +1 snew+t joe +1]. Max jap] Ar (18). 
1 
Wir nehmen nun an, die Abschätzungen (16) seien für alle » bis »= p bewiesen. Dann folgt 
aus (18,) wegen n < c unmittelbar & 
lvr+1| <certier+1, 
d.h. die Richtigkeit für» = p + 1 hinsichtlich v. Dann läßt sich aber eine Schranke für \ 
+1 j 
ty 
1 

angeben: N 
Re CE : 
Berta Ic ER (19), | 
n 1l—ce i 


; 1 x 
wobei ces 5 gemäß (15) und (V) berücksichtigt ist. Weiter gilt aber auch 


1 


‚ c 
2ce S2neNax(1l,n)S; 3 en 
q 2Max (1,2/9) 4° 
so daß für (18,) nunmehr 
3 
z 1w? rl snegetilı (Max IxpI).ep #1 ons a Te er (20) 
geschrieben werden kann. Da analog zu (19) 
2 a p 
Max |oP| = Max (! u) set sich are (21) 
l I 1 1 a 


gilt, ist die rechte Seite in (20) nicht größer als 


n 
pP +1 an r 
CP TI BR (2 +28), 
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so daß wegen 


n 1 1 3g 
— + 28c S——— al, 
ET AE Sm " ZMER(L2f) A 


wPr +1 — Max jwP +1] <eotlgp+l 
ı 


auch 


richtig ist. Damit ist der aufgestellte Hilfssatz bewiesen, denn die suczessiven Approximationen 
konvergieren, und die Lösung liegt in dem Gebiet 


lol, = See Dec 
1 


wie behauptet war. 


3. 
In dem durch die Ungleichungen (13) festgelegten Gebiet existiert aber auch nur eine 
Lösung des Systems (10), wie sich aus dem nächsten Hilfssatz ergibt: 


Hilfssatz 2: Unter den Voraussetzungen des Hilfssatzes 1 existiert in dem durch die Un- 
gleichungen (13) definierten Gebiet nur eine Lösung der Gleichungen (10). 


Wären nämlich ol, a7} und 02,7 zwei verschiedene Lösungen, für welche (13) besteht, so 
müßte die Differenz 


Be En 
= irEenPr 
(A —A + eo) Bı = 3% Mn er Pr — Pic 


genügen. Wegen |o!| < - q finden wir analog zu (18) 


den Gleichungen 


4 || zen Max |ß}| + lo| Max |of| 


und bei Berücksichtigung vonen = & q jedenfalls 


1 
71 Max |ßı] < Jo Max [at] 
und unter Heranziehung von (13) für of: 

1 2 
Max [Bi] S lol 8 ne Max ( n) Be oa 
1 


Andererseits folgt aus (22,) mit (12) 
lol & n e Max |ß}| 
und somit 


2 
lo] = lo] - 8 n2 e2 Max (1, Te 
q q 


auf der rechten Seite durch 1/4 majorisieren läßt, erhalten wir in 


Da sich der Faktor von |o 
jo| = z |o| einen Widerspruch für o # 0. Beio = 0 folgt aber wegen (23) auch , = 0, womit der 
Satz bewiesen ist. 

4. 


Wir wenden uns nun dem Beweis der in der Einleitung gemachten Aussagen zu. Grund- 
lage ist das Bestehen der Voraussetzung (V). Mit Hilfsatz 1 haben wir gezeigt, daß um jeden 
Näherungswert A; ein Kreis vom Radius r =2neMax (1, 2/g) existiert, in welchem mindestens 


ein Eigenwert von A liegen muß. Da die Punkte A,, 3, . . .,A„ in der komplexen Ebene einen Min- 
destabstand g besitzen, der Radius der Kreise wegen (V) aber 


2Max (1,2/Q) 4 


erfüllt, sind diese Kreise punktfremd, d.h. die Eigenwerte von A sämtlich einfach. Damit ist 


der erste Teil der Aussage bewiesen. 
DDE) 


4 


IH 
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= 
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Pr 


H.-G. GISPERT, Über die Nickschwingungen von Flugkörpern bei verschiedenen Geschwindigkeitsgesetzen 395 
mp my  mmmmmthmm  Tn 


ZAMM 41 (1961) Heft 7/8, Seite 3235—330 


Über die Nickschwingungen von Flugkörpern 
bei verschiedenen Geschwindigkeitsgesetzen 


Von H.-G. GISPERT 


Es wird das Schwingungsverhalten von Flugkörpern mit einem Freiheitsgrad bei verschiedenen 
Geschwindigkeitsgesetzen behandelt. Für eine Bewegung mit positiver Beschleunigung erhalten wir 
Amplitudenabnahme. Dagegen können im Fall einer konstanten negativen Beschleunigung (Verzögerung) 
die Schwingungsamplituden nach anfänglicher Abnahme wieder anwachsen. Mathematisch werden diese 
Schwingungsvorgänge i. a. durch höher transzendente Differentialgleichungen charakterisiert. 


The oscillations of a projectile with one degree of freedom are treated. For a motion with positive 
acceleration we get a decrease of amplitudes. However in the case of deceleration the amplitudes of the 
oscillation can increase, when there was a decrease in the beginning. In the mathematical sence these 
oscillations are generally characterized by higher transcendental differential equations. 


PaccmaTpuBaetca KosnebaTelbHbIl IPOoNece MNeTATeIIbHBIX AmmapaToB C OMHOÜ CTEIIEHBIO 
CBOÖOAEI IPH Pa3HbIX 3AKOHAX CKOpocTeü. Ilpım MBu;ReHuMm C IIOJIOFKUTEJIBHEIM YCKOPEeHNEM 
moJIyyaeTrcH yÖbIBaHHe ammauTya. B ciy4yae NOCTOAHHOTO OTPHNATEIBHOTO YCKOPeHHA AM- 
INIAHTYAbI KOJeÖAHHA MOTYT HOCHEe HAYAJIBHOTO YObIBAHHUA CHOBA BO3pacTaTb. Maremarnyecku 
3TU IPOoNecchI xapakrTepusylotca MMdhhepeHmmasbHhIMN YpaBHeHUAMU BbIcleii TPaAHCHeH- 
MEeHTHOCTM. 


1. Einleitung 


Wir denken uns einen Flugkörper so in einer durch seinen Schwerpunkt gehenden Quer- 
achse gelagert, daß die Flugkörperlängsachse in der Vertikalebene (s. Bild 1) um den Schwer- 
punkt drehbar ist. Wird nun der Flugkörper durch einen Luftstrom angeblasen, dessen Richtung 
nicht parallel zur Flugkörperlängsachse ist, oder wird der Flugkörperschwerpunkt geradlinig 
durch die ruhende Luft bewegt, wobei am Anfang der Bewegung die Flugkörperlängsachse um 
einen kleinen Winkel gegenüber der Schwerpunktsbahn verdreht ist, so beginnt der Flugkörper 
um seinen Schwerpunkt zu schwingen und die Amplituden können wachsen, abnehmen oder gleich 
bleiben. 


S: Schwerpunkt, D,: Druckpunkt des Hecks, 
D,: Druckpunkt der Spitze, D,:besamtdruckpunkt 


Bild 1 


Betrachten wir nun den besonderen Fall, daß die Geschwindigkeit des Luftstroms oder 
die des Flugkörperschwerpunktes gleichförmig ist, so lassen sich die Schwingungen mathematisch 
mit Hilfe einer linearen Differentialgleichung zweiter Ordnung mit konstanten Koeffizienten 
beschreiben. Liegt der Druckpunkt des Flugkörpers zudem hinter dem Schwerpunkt, eine Eigen- 
schaft, die hier durchweg vorausgesetzt wird und die man z. B. durch eine entsprechende Flügel- 
fläche am Heck des Flugkörpers realisieren kann, so nehmen die Schwingungsamplituden ab, 
wir erhalten also einen dynamisch stabilen Flugkörper. 

Nehmen wir jetzt einmal an, uns seien die Luftkräfte bekannt, die an einem Flugkörper 
angreifen, dessen Schwerpunkt geradlinig beschleunigt bewegt wird, weiterhin sei es gestattet, 
diese Kräfte mit Hilfe des Newroxschen Ansatzes darzustellen, so lassen sich die durch eine An- 
fangsstörung verursachten Schwingungen der Körperlängsachse um die Richtung der Schwer- 
punktsbahn auch wieder durch eine lineare Differentialgleichung zweiter Ordnung charakterisie- 
ren. Im Gegensatz zur oben erwähnten Differentialgleichung sind die Koeffizienten dieser Glei- 
chung nicht mehr Konstante, sondern Funktionen der unabhängigen Veränderlichen. Als 
Lösungen ergeben sich i. a. höher transzendente Funktionen mit oszillierendem Charakter. Unter 
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Berücksichtigung der bereits erwähnten Einschränkungen werden in der vorliegenden Arbeit 
Schwingungen eines Flugkörpers untersucht, dessen Schwerpunkt nach dem Geschwindigkeits- 
gesetz A B 
v=—kt—k, 

geradlinig bewegt wird. Darin ist k, eine positive und k, eine beliebige Konstante. Für eine Bewe- 
gung mit positiver Beschleunigung erhalten wir Amplitudenabnahme. Dagegen können im Fall 
einer konstanten negativen Beschleunigung (Verzögerung) die Schwingungsamplituden nach 
anfänglicher Abnahme wieder anwachsen. 


2. Momentengleichung 


Wir legen den Koordinatenursprung in den Flugkörperschwerpunkt. Weiterhin nehmen wir 
an, daß der Anstellwinkel « auf kleine Beträge beschränkt bleibt, damit wir, aus Gründen der 
Vereinfachung, einen linearen Ansatz erhalten. Für konstantes Trägheitsmoment I ergibt sich 
dann die Momentengleichung, welche die Nickschwingungen der Flugkörperlängsachse charak- 
terisiert, (s. Bild 1) u 

IA=—N,.L+HNL #2 22 2, Due 


mit N,: Normalkraft, die auf den beflügelten Teil des Flugkörpers wirkt, 
N,: Normalkraft, die auf die Flugkörperspitze wirkt, 

l;: (i = 1,2): die entsprechenden Druckpunktabstände. 

Führen wir noch die Voraussetzung ein, daß die Änderung von N, bzw. N, in bezug auf 
den Anstellwinkel, entsprechend dem Newroxschen Ansatz, proportional zum Staudruck ist, so 
können wir schreiben 

N=-htwrla 4 53 
ee SR, m: 


2, ) 


Darin sind cy, = konst. die entsprechenden Normalkraftbeiwertsänderungen, og die Luftdichte 
und F eine gemeinsame Bezugsfläche. Mit diesen Ansätzen geht (1) unter der Bedingung von 


Ne ner (=) 


h eh ch, L— ceul 
in ; 
: + i 5 ® 
I&= —{cy, +} 5 Foa—cyh5 Fura ae Se 
über; 
mit h: Abstand des Gesamtdruckpunktes vom Schwerpunkt, 
cx: Normalkraftbeiwertsänderung des Flugkörpers. 
Setzen wir nun 


I MET ORE® ‚ 
at nMSF: bh- 
so ergibt sich aus (2) 


Star bRamÜ,. una er 


Sind a,, b, und v Konstante, so setzt sich bekanntlich die allgemeine Lösung von (3) aus 


einfach transzendenten Funktionen zusammen. Für die nun folgenden Untersuchungen nehmen 
wir zunächst an 


a,, b,: konst. ; v: nicht konst. 
Daraus ergibt sich, daß schon für einfache Geschwindigkeitsgesetze die Lösungen von (3) 


aus der Klasse der einfach transzendenten Funktionen in die der höheren, z.B. der BESSEL-, 
der konfluenten hypergeometrischen und der hypergeometrischen Funktionen übergehen. 


3. Aufstellung der Geschwindigkeitsgesetze 


ee an Bewegungsgleichung des Schwerpunktes eines Flugkörpers in Richtung der Flugbahn 
autet 


mim=S—W-—Gr, 
wobei m die Masse des Flugkörpers, S der Schub, W der Luftwiderstand und Gr die Komponente 
des Gewichtes in Flugbahnrichtung ist. Setzen wir nun 


Ww S nn G 
ke und a Kon k,: konst. , 


4 
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so leiten sich aus der Beziehung 
die hier betrachteten Geschwindigkeitsgesetze ab. 
Die Integration von (4) liefert 


= ts jetz En Var. EEE 6), 


mit v(0) = v,. Für k,— 0, v, = 0 und k, < 0 ergibt sich 
v=—kt, A U N & (6), 
d. h. der Flugkörper erhält eine konstante Beschleunigung, die etwa mit Hilfe eines Raketentrieb- 


werkes erzeugt wird. 
Für k, —0 und k, > 0 geht (5) in 


v=-krt, ee a ee EN ER (7) 
k, 


über, d. h. der Flugkörper erhält eine konstante Verzögerung, die etwa durch Bremsraketen ver- 
ursacht wird. 


Für k, —0 und k, > 0 folgt schließlich aus (5) 


d. h. der Flugkörper erhält eine Verzögerung, die sich z. B. infolge eines Luftwiderstandes ergibt. 


4. Beschleunigte Bewegung 
Mit 6, a=—qak,b=b,k% geht (3) in 


OPER OA DR ae DE ee (9) 
über. 


In dem Spezialfall a = 0 [1] erhält man mittels der Variablen-Transformation z = in R 


[72 
und der Substitution x = /z- u die Bzsseusche Differentialgleichung 


du du 1 
Zee Them 


deren allgemeine Lösung sich aus einer Linearkombination der Bessefunktionen mit den Zeigern 
v = 1/4 und v = — 1/4 zusammensetzt. 


Fordern wir als Anfangsbedingung «(0) # 0, &(0) = 0, so ist 


= C- Vz J_1/4(2) . 
Betrachten wir nun den allgemeinen Fall, so geht (9) mit Hilfe des Produktansatzes 


= en®.yu, 
worin x, eine Lösung von 
TE EA (10) 
ist, und mit 


A=4m,+a, Bit, Ee= — — 


deu 1 du = ‚B 


über. Das ist die konfluente hypergeometrische Differentialgleichung [2], [3] mit der allgemeinen 
Lösung E E 
Da! 1 
u kr Di RT 1/2 BI . 
= Aaa)+e 2 le ea) 
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Unter Berücksichtigung der Anfangsbedingungen a(0) # 0, &(0) = 0 und 
ee 
N Pre er dere) 


erhalten wir mithin 
227, ta w 


2x, ta 1 4. +a ) 
Bee ESF NN = ec) Re 1... 11). 
Fler 2 aD 


Wenn man andererseits versucht, (9) mit Hilfe des Reihenansatzes 


Velne . 


s= a, 
nu=0 
zu lösen, so führt das auf die Rekursionsformel 
ir er. (12). 
a_,„=0 


Aus (11) erkennt man, daß «(f) das Produkt von zwei ganzen Funktionen ist, mithin a(f) 
selbst eine ganze Funktion ist. Zur geometrischen Deutung von «(f) führen wir in (9) den Ansatz 


uote 0 ee a (13) 
ein und erhalten 
ü+ ER Be ne; 
ee 


Da im allgemeinen für Flugkörper a < b ist, können wir für hinreichend große t 
2 
ü+ (bg) Fu =0 


setzen. Die Lösung dieser Differentialgleichung kann wieder, bei Vorgabe entsprechender An- 
fangsbedingungen, auf die BesseLfunktion mit dem Zeiger — 1/4 zurückgeführt werden. Mithin 


erhält man 
1 a2 
u= (* VEI_i = \; euT P) 


.(C*: konst., von a, b abhängig) 
und für hinreichend große f? ergibt sich die asymptotische Lösung [4] 
1 11/5 @ e) 
D.C#* — en —— P——) 2.2... 0... ; 
us.G 7 cos(4 /5 7! 3 et. 
(C**: konstant; von a, b abhängig). 


Wir entnehmen daraus folgendes: Die Lösungsfunktion u oszilliert. Durchlaufen wir die 


positive Zeitachse i von links nach rechts, so nehmen die Amplituden wie 1//t ab. Wegen £ im 
Argument der cos-Funktion verkleinern sich mit wachsendem ? die Nullpunktabstände. Daraus 


a 


| 


folgt, daß wegen (13) die Amplituden von «x(!) wie e * 5 abnehmen. 
ER 


Führen wir schließlich in (9) als neue Veränderliche den Flugweg s (englisch: traveled distance) 
ein, so ergibt sich 
dx 


C at 
retratee-t 


€ € Cs konstant. Für hinreichend große I erhalten wir somit die Schwingungsgleichung mit kon- 
stanten Koeffizienten. 
5. Verzögerte Bewegung 
Mit (7) und wegen a(l) = x(— I) ergibt sich jetzt aus (3) 
s— ats +bRa=0. 


Diese Gleichung!) ist auch vom Typus der konfluenten h ri i l 
j ypergeometrischen Differential- 
gleichung. In (10), (11), (12), (13) und (14) ist also a durch — a zu ersetzen. Aus (14) erkennen 


!) Es ist wieder 7 = t gesetzt. 


BE Nenn 


nn u en BEE 
N je 
4 ers u ' De: . 


rin WEN 
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abnahme, die entweder für t < !* fortgesetzt wird oder aber 


m — 
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wir wieder, daß die Amplituden von u(!) mit wachsendem t wie 1//t abnehmen, dagegen nehmen die 
von «(f) wegen des divergenzerzeugenden Faktors exp (a 12/4) von einer bestimmten Stelle {* 
an mit wachsendem ? zu. Für t< {* kann es ein endliches Zeitintervall geben, in dem die Ab- 
nahme der Amplituden von u(f) stärker ist als das Wachstum 
von exp (a2/4). | 

Mithin erhalten wir folgenden Schwingungsverlauf: 
Wegen (7) und a(t) = a(— f) müssen wir die positive l-Achse 
von rechts nach links durchlaufen. Es ergibt sich zunächst bei 
abnehmendem { bis zum Zeitpunkt {* eine Amplituden- 


in eine Amplitudenzunahme übergeht, obwohl wir hier nur 
solche Flugkörper zugelassen haben, die im stationären Flug Bild 2 


dynamisch stabil sind. Bild 2 skizziert den Schwingungsverlauf. 


Mit (8), d. h., wenn der Flugkörper nur durch den Luftwiderstand ‚verzögert wird, 
t=uwkt+l, K,=alk, K,=b,/k geht (3) in die EuLersche Differentialgleichung 


BR A FR 
dee Mae 


-I- 


über. Die Lösung lautet 


ge Aa WER ar 
edel iu, kt +1} ‚co —( 5 . MONee): 


Be, wer: 
+6w,kt+1} ? »sin Ir-( 5 ) In (w,kt-+ ” air 
Daraus erkennen wir: Für wachsende Z, {> 0, nehmen die Amplituden von a(f) entweder 


nur zu (0 < (1 — K,)/2) oder nur ab (( > (1 — K,)/2). 
Führen wir nun (5) in (3) ein und substituieren 


1-+tg {ar tg Vr v—Vkı-ko- N = + ee ri (16), 
2 


so erhalten wir die hypergeometrische Differentialgleichung 


Sa u 
mit der allgemeinen Lösung 
1 
— [kk—4,+4,] 1 1 1 
a CE P (gg lat + Ali gg hr AA gpA+lid) 
aa (1 1 1 
ee EB mare äl; Tee: on &) (17); 


A, = VI — a” —4Ab,; A=V®—4Ab,. 


Für k,— 0 erhält man, wenn in (17) C, und (C, geeignet gewählt werden, wieder (15). Ebenso 
ergibt sich für k, —0, falls man & = 1 — n setzt, wieder (11). Für entsprechende k, und k, kann 


man mithin die Schwingungsformen von (17) charakterisieren. 


6. Eine Verallgemeinerung der Ausgangsgleichung 


Nehmen wir jetzt an, daß alle Parameter, die in den Koeffizienten von (3) auftreten, mit 
der Zeit £ veränderlich sind, z.B. I = /{t), I, = L(f) usw., so erhält die Momentengleichung die 
allgemeine Gestalt 

ee Le (18). 
Mit 


ergibt sich 


&. KK u EZ Ze j 
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Der folgende Oszillationssatz [5] gestattet uns, Aussagen über u(f) zu machen. 


In 
y’ + go) y) = 0 


sei p(x) stetig und monoton wachsend für x > x, /(y) stetig, monoton wachsend, eine ungerade 
Funktion und erfülle für |y| < Y, eine Lipscurrz-Bedingung. Ist y, mit 0 < |y,| < Y, gegeben, 
so existiert die durch die Anfangswerte y(X,) = Yo Y’(X%) = 0 bestimmte Lösung für alle x > x, 
hat unendlich viele Nullstellen und die Amplituden nehmen monoton ab, haben aber nicht not- 
wendig den Limes 0. 

Ein Vergleich mit (19) liefert, daß u die für /(y) geforderten Bedingungen erfüllt. Genügt 
nun andererseits (4 n — m? — 2 m)/4 den Voraussetzungen von p(x) und ist außerdem - 


t 
lim [mdt= co, 
t»+ohb 
so gilt für die Lösung von (18) 
lim oa) =0, E 


t+ 00 


d. h. die Amplituden konvergieren nach Null. 
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L. Collatz, The Numerical Treatment of 
Differential Equations. XV + 5688. m. 118 
Abb. Berlin/Göttingen/Heidelberg 1959. Springer- 
Verlag. Preis geb. 93,60 DM. ü 


Die vorliegende englische Ausgabe des bekannten 
Standardwerkes über die numerische Behandlung 
von Differentialgleichungen ist im wesentlichen eine 
Übersetzung der zweiten deutschsprachigen Auflage, 
die in dieser Zeitschrift 35 (1955), S. 320, ausführlich 
besprochen worden ist (Besprechung der ersten Auf- 
lage: ZAMM 31 (1951), S. 395). Selbstverständlich 
sind die notwendig oder zweckmäßig erscheinenden 
Verbesserungen und Ergänzungen vorgenommen wor- 
den. Vor allem sind den Literaturhinweisen neue 
hinzugefügt und die an sich schon reichhaltige Samm- 
lung durchgeführter Beispiele ist durch neue Auf- 
gaben erweitert worden. Das bewährte Werk ist 
nicht nur wegen der Darstellung der Methoden, 
sondern vor allem wegen der gründlichen Behandlung 
der Fragen der Fehlerabschätzung längst zum zu- 
verlässigen und unentbehrlichen Ratgeber aller der- 
jenigen geworden, die Differentialgleichungen nume- 


‚ risch zu behandeln haben. Es ist zu hoffen, daß die 


Übersetzung ins Englische nur dazu dient, dem Buch 
die verdiente Verbreitung in das englische Sprach- 
gebiet zu verschaffen, daß sie aber nicht die deutsche 
Ausgabe ersetzen soll. 


Dresden H. HEINRICH 


D. Kappos, Strukturtheorie der Wahrschein- 
lichkeitsfelder und -räume. IV + 1368. Berlin/ 
Göttingen/Heidelberg 1960. Springer-Verlag. Preis 
brosch. DM 21,80 DM. 


Aus logischen und formalen Gründen erscheint es 
angemessen, bei Vorliegen eines bestimmten Ver- 
teilungsgesetzes p auf einer o-Mengenalgebra die 
Ereignisse der Wahrscheinlichkeit Null mit dem un- 
möglichen Ereignis zu identifizieren, d.h. zur Rest- 
klassenalgebra B nach dem o-Ideal der Ereignisse der 
Wahrscheinlichkeit Null überzugehen. Da p(A) nur 
von der Restklasse Ä von A abhängt, dürfen wir 
p(A) = q(Ä) setzen und erhalten so ein o-Wahr- 
scheinlichkeitsfeld [B, q], d.h. ein geordnetes Paar, 
bestehend aus einer o-vollständigen Booteschen Al- 
gebra und einem strikt positiven Verteilunggesetz über 
dieser Algebra. Die Darstellung von Karros beginnt 
mit dem schwächeren Begriff des Wahrscheinlich- 
keitsfeldes, in dem auf die o-Vollständigeit von B so- 
wie die o-Additivität von g verzichtet wird. Durch 
Adjunktion ‚idealer‘ Ereignisse kann man stets zu 
einem o-Wahrscheinlichkeitsfeld übergehen, das durch 
einfache Forderungen bis auf Isomorphie eindeutig 
bestimmt wird. 

Zur Strukturanalyse von o-Wahrscheinlichkeits- 
feldern wird der Begriff der stochastischen Unab- 
hängigkeit von Familien von o-Unteralgebren benutzt, 
wobei sich wegen der strikten Positivität der Wahr- 
scheinlichkeit besonders einfache Zusammenhänge 
ergeben. Die Darstellung führt bis zur MAHnARAMschen 
Klassifikation. 

Besonders interessant ist Kapitel VII, in dem die 
Theorie der kompakten Maße und ihre Anwendung 
auf die Konstruktion von Maßen in Produkträumen 
dargelegt wird. Das Buch endet mit einem Bericht 
über das von REnYt eingeführte, auf dem Begriff der 
bedingten Wahrscheinlichkeit beruhende Axiomen- 
system. 

Im Vorwort spielt der Verfasser darauf an, daß der 
„punktfreie‘“ Aufbau der Wahrscheinlichkeitsrech- 
nung bestimmte technische Vereinfachungen — z. B. 
die Aufhebung von Mehrdeutigkeiten — mit sich 


Buchbesprechungen 
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bringt. Leider werden diese Fragen jedoch nicht dar- 
gestellt, wie überhaupt der Zusammenhang zur Wahr- 
scheinlichkeitsrechnung recht locker ist. 


Strausberg b. Berlin K. MATrTHEs 


H.-U. Smoltezyk, Ermittlung eingeschränkt 
plastischer Verformungen im Sand unter 
Flachfundamenten. 1108. m. 52 Abb. u. 4 Tab. 
Berlin 1960. Verlag Wilhelm Ernst & Sohn. Preis 
brosch. DM 16,80. 


Der erste Teil der Arbeit beschäftigt sich mit 
theoretischen Untersuchungen über den Grundbruch 
in einem homogenen isotropen Boden. Ausgehend 
von einem einfachen und anschaulich abgebildeten 
Scherversuch wird der Spannungszustand im Boden 
unter Voraussetzung der CovLomgschen Bruchtheorie 
für den als ideal-plastisch angenommenen Sand defi- 
niert. Die Gleitlinien werden aus einzelnen Kurven- 
stücken zusammengesetzt. Betrachtungen eines verti- 
kal bzw. schief belasteten Flachfundamentes auf 
„gewichtslosem‘‘ Boden liefern die mathematische 
Form der Gleitlinien bzw. die Randbelastung. Die 
ausführlichen Ableitungen der entsprechenden Be- 
dingungen für ‚schweren‘ Boden werden durch je 
ein Zahlenbeispiel und die Abbildung eines Labor- 
versuches erläutert. Ein Vergleich der gewonnenen 
Ergebnisse mit der herkömmlichen Grundbruchberech- 
nung und praktischen Erfahrungen zeigt keine be- 
friedigende Übereinstimmung, was auf Überlagerung 
kleinerer Brucherscheinungen, Begrenzung der freien 
Beweglichkeit der Fundamente und die meist vor- 
handene räumliche Tragwirkung zurückgeführt wird. 

Im zweiten Teil werden eingeschränkt plastische 
Verformungen von Sanden untersucht, worunter die 
Deformation infolge Verdrängung und Verdichtung 
bzw. Auflockerung verstanden werden. Die aus der 
dabei entwickelten Theorie an Hand eines Zahlen- 
beispieles gewonnenen Ergebnisse werden mit voran- 


‚gestellten Versuchsergebnissen verglichen und zeigen 


gute Übereinstimmung. 

Für die Praxis besitzen die beiden Verfahren in- 
folge ihres großen mathematischen und rechneri- 
schen Aufwandes noch keine Bedeutung. Die Arbeit 
ist jedoch allen zu empfehlen, die sich mit Forschun- 
gen über Spannungsverhältnisse, sowie Bruch- bzw. 
Fließbedingungen in Baugrund befassen, zumal nicht 
nur die mathematischen Ansätze ausführlich be- 
handelt werden, sondern auch aus den Ergebnissen 
üblicher Laborversuche recht interessante Schlüsse 
gezogen werden. 


Dresden H. E. KLAHN 


C. Jacob, Introduction Math6matique & la 
Mecanique des Fluides. 1286 S. m. 258 Abb. 
Paris 1959. Gauthier-Villars. Preis geb. $ 10,—. 


In diesem umfangreichen Werk, das 1952 in rumä- 
nischer Sprache erschien und dessen französische 
Übersetzung und Neubearbeitung jetzt vorliegt, 
unternimmt der Verfasser den Versuch, drei Gebiete, 
die im allgemeinen in getrennten Lehrbüchern dar- 
gestellt werden, möglichst vollständig unter einheit- 
lichen Gesichtspunkten in einem Band zu behandeln. 
Es handelt sich dabei erstens um die mathematische 
Theorie der Lösung von Randwertproblemen 1. und 
2. Art für die Potentialgleichung mittels konformer 
Abbildung und durch Integralgleichungen, zweitens 
die Theorie der inkompressiblen und drittens der 
kompressiblen Strömungen. Abgesehen von der Her- 
leitung der NAvIER-Stokzsschen Gleichungen, der 
Erwähnung der Stokzsschen Näherung und der Dar- 
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stellung einiger einfacher Lösungen, wurden nur 


zähigkeitsfreie Strömungen behandelt. Das Wort 


„Introduetion‘‘ im Titel des Buches ist irreführend, 
denn es ist als einführendes Lehrbuch ungeeignet, da 
der Verfasser auf jede physikalische Diskussion der 
Voraussetzungen und Ergebnisse und die Berechnung 
von Beispielen verzichtet. Ein unerfahrener Leser 
bekommt also überhaupt kein Gefühl für den Gültig- 
keitsbereich der verschiedenen Lösungen und für die 
Vor- und Nachteile der einzelnen Näherungsansätze. 
Das Buch ist dagegen als Nachschlagewerk für den 
Wissenschaftler sehr zu empfehlen, denn in zahlreichen 
Gebieten werden auch noch die neuesten Ergebnisse 
(teilweise bis zum Jahre 1957) berücksichtigt. Es ist 
schwer, eine genaue Inhaltsangabe zu machen. Das 
Buch enthält wohl alle wichtigen Grundlagen über 
konforme Abbildung, Potentialtheorie, Profil- und 
Tragflügeltheorie, Wirbelsätze, Singularitäten, Stoß- 
wellen, Hodographenmethoden und Näherungsver- 
fahren für die Berechnung von Unter- und Überschall- 
strömungen, aber die Auswahl und der Umfang in der 
Behandlung der speziellen Probleme ist naturgemäß 
subjektiv durch die Interessen des Verfassers bestimmt. 
Dem Referenten fiel etwa auf, daß bei der Behandlung 
schallnaher Strömungen noch die neueren Ergebnisse 
und die Kontroverse über lokale Überschallgebiete 
erwähnt werden, und daß das Buch einen aus dem 
übrigen Zusammenhang etwas herausfallenden Ab- 
schnitt über das wenig behandelte Gebiet der Scher- 
strömung hat, während andererseits in dem Kapitel 
über Freistrahlen und Totwassergebiete jeder Hinweis 
auf die wichtigen neuen’ Arbeiten über endliche Tot- 
wasser mit Unterdruck fehlt. Oberflächenwellen 
werden gar nicht behandelt, Gitterströmungen nur 
sehr kurz. Sehr wertvoll sind die zahlreichen Litera- 
turhinweise, insbesondere werden auch viele Arbeiten 
aus Rumänien und der Sowjetunion zitiert, die noch 
vielfach unbekannt sein dürften. Das Buch hat leider 
kein Sach- und Namensregister, aber nach kurzer Ein- 
arbeitung gestattet das umfangreiche Inhaltsverzeich- 
nis ein rasches Auffinden jeder gesuchten Stelle. 


Freiburg i. Br. G. JUNGCLAUS 


Combinatorial Analysis. (Proceedings of the 
Symposia in Applied Mathematics, Vol. 10) VI+ 
311 S. Providence, Rhode Island, 1960. American 
Mathematical Society. Preis geb. $ 7,70. 


Aus der Zahl der 24, im vorliegenden Bande ver- 
einigten Artikel mögen hier nur einige herausgegriffen 
werden, bei denen die Beziehung zur Numerik bzw. 
zur Anwendung besonders ins Auge fällt; ein Wert- 
urteil über die nicht erwähnten Titel soll damit keines- 
wegs ausgesprochen werden. 


R. BELLMAN, der zusammen mit H. Haus JR. als 
Herausgeber dieser Artikelsammlung zeichnet, stellt 
unter der Überschrift ‚„Kombinatorische Prozesse 
und dynamisches Programmieren“ seine Methode der 
rekursiven Funktionalgleichungen zur Lösung be- 
stimmter Typen von Programmierungsaufgaben dar. 
Diese Methode ist sowohl im Falle kontinuierlicher 
wie auch diskreter(!) Variablen anwendbar. Er skiz- 
ziert, über sein Buch ‚Dynamic Programming‘ hin- 
ausweisend, numerische Kunstgriffe zur effektiven 
Berechnung der rekursiv definierten Funktionen und 
kündigt weitere Veröffentlichungen an. Viele kombi- 
natorische Optimierungsprobleme lassen sich auf 
„ganzzahlige‘ lineare Programmierungsprobleme zu- 
rückführen ; für solche Probleme beschreibt R. GoMoRY 
ein Verfahren, welches endlich viele Male die Anwen- 
dung der Simplexmethode verlangt. Der Artikel von 
R. Karapi behandelt Verteilungsaufgaben in Fern- 
sprechnetzwerken. M. GERSTENMAYER gibt Rat- 
schläge zur Lösung von Transportproblemen großen 
Umfangs. M. FLooD präsentiert einen neuen Beweis 
für den bekannten Satz von Könıg-EGERVÄRY und 
gründet darauf eine Modifikation der „ungarischen 


Methode“ 


"angestellt, wie man einen Rechena maten { 


> 


ee 


problems (ganzzahliges Traı oble 
b; =»). Te Tausche koche Verall; ın 
sung des Transportproblems mit. 21,421 


cheint damals noch nicht bekann 
schein Reihenfolge 


J. FOULKES ie den 

Einteil er Knotenpi  zugeor 
phen ne Äquivalenzklassen beizukomm« 
MANN gibt weitere Beispiele für die erfolgr 
nutzung der Theorie des linearen Pro nmierens 
der Klärung mengentheoretischer Fragen. 
Artikeln von WALKER, PAIGE, TOMPKINS, GLEAS 
LEHMER und Swırr werden wertvolle Üborieen 


io. 
Lösung kombinatorischer Fragen einsetzen kann. Als 

Autoren sind noch zu nennen: KuHN, Tucker, Hart, 
Bruck, HuGHEs, ALBERT, DILWORTH, WHITEMAN, 
Ryser, Marcus and NEwMAN, TAUSSKY, Topp und 
GILBERT. Ein Teil der Arbeiten befaßt sich allerdings 

auch mit Problemen der „abzählenden‘‘ Geometrie, 

der Algebra und Zahlentheorie. 

Die wärmstens empfohlene Lektüre des vorliegen- 
den Sammelbandes wird für viele angewandte Mathe- 
matiker ein Tor in ein unbekanntes, aber sehr reizvolles 
und z. T. unerforschtes Land aufstoßen. 


Dresden L. BITTNER 


” 


W. Flügge, Stresses in Shells. XII + 4998. 
m. 244 Abb. Berlin/Göttingen/Heidelberg 1960. 
Springer-Verlag. Preis geb. DM 58.80. 3 


Bereits von 25 Jahren hatte der Verfasser ein 
Buch über Schalentheoriein deutscher Sprache heraus- 
gegeben. Durch seine Lehrtätigkeit nach dem zweiten 
Weltkrieg an der Stanford-University in Californien 
als einer der Nachfolger TIMosHENKoS reifte in- 
zwischen sein Entschluß, ein Schalenbuch in engli- 
scher Sprache herauszugeben, welches nunmehr vor- 
liegt. Zwischen den beiden extremen Möglichkeiten, 
entweder die mathematische Seite der Schalentheorie 
in den Vordergrund zu stellen und praktisch ein 
Buch über die Behandlung spezieller Differential- 
gleichungen zu schreiben oder den Inhalt ganz der 
Anwendung auf Schalendächer und Druckgefäße mit 
Berücksichtigung aller rechnerischen und konstruk- 
tiven Einzelheiten zu widmen, hat Verfasser offenbar 
einen glücklichen Mittelweg eingeschlagen. Die 
mathematischen Voraussetzungen entsprechen ganz 
dem jeweils behandelten Problem; teils genügen 
Differentialgleichungen mit konstanten Koeffizienten, 
teils solche mit veränderlichen Koeffizienten, Produkt- 
lösungen oder komplexe Funktionen. 

Das Buch zerfällt in 4 Hauptabschnitte. Der erste 
Abschnitt (Kap. 1) liefert die grundlegenden Begriffe 
und Definitionen. Der zweite Hauptabschnitt (Kap. 
2—4) ist der Theorie der Membranschalen gewidmet. 
Im Hinblick auf die praktischen Erfordernisse sind 
verschiedenartige Schalenformen und Auflagerbedin- 
gungen behandelt. Durch viele Beispiele wird der 
Rechnungsgang bei Festigkeitsproblemen für Flüssig- 
keitsbehälter und Schalendächer erläutert. Kritische 
Betrachtungen zeigen den Gültigkeitsbereich der 
Membrantheorie. Der dritte Hauptabschnitt (Kap. 5 
u. 6) behandelt die Biegetheorie der beiden wichtig- 
sten Schalentypen: Kreiszylinderschale und Rota- 
tionsschale. Verf. hat sich hierbei bemüht, die mathe- 
matische Behandlung so übersichtlich und einfach 
durchzuführen, wie es mit Rücksicht auf die Genauig- 
keit der Theorie zulässig ist. Deshalb wurde auf die 
Anwendung tensorieller differentialgeometrischer Me- 
thoden, wie sie sich in der allgemeinen Schalentheorie 
eingebürgert haben, verzichtet. Für eine große Anzahl 
der behandelten Probleme ist die Lösung bis zur An- 
gabe direkter Anwendungsformeln vorgetrieben, wel- 
che die wesentlichen Kraft- und Verschiebungsgrößen 
liefern. Vielfach sind die erforderlichen Parameter 
in Tabellenform angegeben. Der letzte Hauptab- 
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schnitt (Kap. 7) behandelt die Stabilitätstheorie der 
Schalen. Die Genauigkeit einzelner Methoden der 
mathematischen Stabilitätstheorie der Schalen ist 
z. Zt. noch teilweise umstritten; zudem zeigt sich oft 
ein verhältnismäßig großer Unterschied zwischen der 
experimentellen Beullast einer Schale und dem ent- 
sprechenden Rechnungswert. In Einzelfällen hat die 
Berücksichtigung endlicher Deformationen eine ge- 
wisse Klärung gebracht; der erforderliche numerische 
Aufwand kann aber dem Konstrukteur nicht zu- 
gemutet werden. Deshalb hat sich Verf. auf eine 
Auswahl solcher Stabilitätsaufgaben beschränkt, wel- 
che z. Zt. als ausreichend geklärt angesehen werden 
können. i 

Das Buch wendet sich in erster Linie an Ingenieure 
der Praxis, aber auch an Studenten der höheren 
Semester. Infolge seiner auf die praktische Anwen- 
dung gerichteten Zielsetzung wird es sicher großen 
Anklang finden. 


München H. NEUBER 


Handbuch der Physik, Band IX. Strömungs- 
mechanik III. Herausg. S. Flügge. VII + 8158. 
m. 248 Abb. Berlin/Göttingen/Heidelberg 1960. Sprin- 
ger-Verlag. Preis geb. DM 198.—. 


Die sechs im vorliegenden Band des Handbüchs 
für Physik vereinten Artikel aus verschiedenen Ge- 
bieten der theoretischen Strömungslehre sind im 
wesentlichen mathematischer Natur, nur gelegentlich 
finden sich Hinweise auf experimentelle Ergebnisse. 
Der Band vermittelt einen Eindruck von dem Aus- 
maß, in dem mathematische Methoden in der Strö- 
mungslehre Eingang gefunden haben. 

Der erste Artikel von M. SCHIFFER handelt von der 
analytischen Theorie der stationären Unter- und 
Überschallströmungen bei Vernachlässigung von Rei- 
bungs- und Wärmeleitungseffekten und (im wesent- 
lichen) Beschränkung auf ideale Gase. Nach der 
Diskussion der Charakteristikentheorie folgen weitere 


Kapitel über die linearisierte Theorie, die Hodo- - 


graphenmethode und die Theorie zweidimensionaler 
Unter- und Transschallströmungen mit Existenz- und 
Eindeutigkeitssätzen. 

Der darauffolgende Artikel von H. CABANNES 
bringt eine Diskussion der in dem ersten Artikel nicht 
behandelten Verdichtungsstöße in stationärer, rei- 
bungsloser Strömung eines idealen Gases, wo Ver- 
dichtungsstöße als Singularitätsflächen auftreten. 
Probleme der Stoßstruktur usw. werden nicht erörtert. 

Ein Artikel von R. E. Meyer über die Charak- 
teristikentheorie der reibungsfreien Gasdynamik be- 
faßt sich mit den analytischen Aspekten dieser 
Theorie und nicht mit den auf ihr aufbauenden 
zahlreichen numerischen Verfahren. 

Es schließt sich ein kurzer Artikel von R. TImMAN 
über die linearisierte Theorie der instationären, rei- 
bungslosen, kompressiblen Strömungen an, der im 
wesentlichen einen Bericht über die instationäre 
Tragflügeltheorie gibt. 

D. GitBarG behandelt in einem Artikel über 
Strahlen und Totwasser die Freistrahltheorie bei 
Vernachlässigung der Reibung und Kompressibilität 
des strömenden Mediums. Dieser spezielle Zweig der 
Potentialtheorie ist von Bedeutung bei der Erklärung 
der wirklich beobachteten Freistrahl- und Totwasser- 
erscheinungen, vor allem aber für die Theorie der 
Kavitation. Nach einem Überblick über spezielle 
Lösungen werden allgemeine Sätze, darunter Existenz 
und Eindeutigkeitssätze mitgeteilt und abschließend 
numerische Methoden erörtert. 

In dem letzten, fast die Hälfte des Bandes füllenden 
Artikel behandeln J. V. WEHAUSEN und E. V. Lar- 
TONE die Theorie der Oberflächenwellen auf einer 
schweren Flüssigkeit unter Beschränkung auf die 
grundlegenden Dinge und Verzicht auf allzu ausführ- 
liche Diskussion von Anwendungen. Den zugrunde- 


liegenden Näherungstheorien entsprechend ist die 
Hauptmasse des Stoffes in einen Abschnitt über 
Wellen infinitesimaler Amplitude und einen Ab 
schnitt über Flachwassertheorie gegliedert. Der erste 
dieser Abschnitte beginnt mit einer Übersicht über 
spezielle Lösungen, besonders solche für zeitlich 
stationäre oder pulsierende quell- oder wirbelförmige 
Singularitäten unter der freien Oberfläche der ruhen- 
den oder bewegten Flüssigkeit. Es folgen: Diskussion 
der Gruppengeschwindigkeit, Wellen infolge eines 


. Hindernisses im Strom, infolge bewegter Begren- 


zungswände, Schiffswellen, Unterwassertragflügel, An- 
fangswertprobleme, Wellen in begrenzten Bassins, 
Einfluß der Oberflächenspannung und der inneren 
Reibung, höhere Näherungen. 

Auf den Abschnitt über Flachwassertheorie, in dem 
neben der linearisierten Theorie die den entsprechen- 
den gasdynamischen Verfahren ganz analoge Charak- 
teristikenmethode der nichtlinearen Theorie sowie 
die Einzelwelle und die KoRTEwEG-DE Vrızs-Wellen 
erörtert werden, folgt ein letzter Abschnitt über 
exakte Lösungen. Dieser enthält außer einigen streng 
gültigen Theoremen, z. B. solchen, die den Massen-, 
Impuls- und Energietransport betreffen, und expli- 
ziten Lösungen, wie z.B. den GERSTNER-Wellen, 
auch einen Überblick über Existenzsätze der Wellen- 
theorie. 

Da die ausführliche Zusammenfassung der Wellen- 
theorie in dem Lehrbuch ‚‚Hydrodynamik“ von 
Lamg (6. Aufl. 1932) nur die ältere Literatur berück- 
sichtigt und da das Buch ‚Water Waves“ von 
STOKER (1957) das Thema stark selektiv behandelt, 
gibt der vorliegende Artikel einen sehr wünschens- 
werten und wertvollen Überblick gerade über die 
neuesten Entwicklungen auf diesem mathematisch 
wie physikalisch gleichermaßen interessanten Gebiet. 


Freiburg i. Br. E. BECKER 


F.R. Güntsch, Einführung in die Program- 
mierung digitaler Rechenautomaten. Mit be- 
sonderer Berücksichtigung der Z 22. 1448. Berlin 
1960. Walter der Gruyter & Co. Preis geb. DM 24,—. 


Die vorliegende Einführung in die Programmie- 
rungstechnik vermittelt hauptsächlich die Beschrei- 
bung der Programmierungseigentümlichkeiten des 
Rechenautomaten Z 22. Der Text ist knapp gehalten, 
dafür sind viele erläuternde Beispiele vorhanden, 
mehr als ein Drittel des Umfanges sind für Fluß- 
diagramme und ausgeführte Programme verwendet. 
Im einzelnen wird zunächst die Beschreibung von 
Rechenprozessen mittels der algorithmischen Schreib- 
weise nach RUTISHAUSER, die Zahlendarstellung in 
elektronischen Rechenautomaten und eine kurze Be- 
schreibung der Z 22 gegeben. Nach Einführung des 
Flußdiagramms folgen einige einfache Program- 
mierungsbeispiele. Danach werden die grundlegenden 
Begriffe der Programmierung, Struktur von Program- 
men, Unterprogrammtechnik und Adressenänderun- 
gen erörtert und jeweils durch Beispiele ergänzt. 
Sämtliche Programme sind in dem für die Z 22 üb- 
lichen Externcode geschrieben, es werden deshalb für 
symbolische Adressen auch nicht beliebige Symbole, 
sondern nur eingeklammerte natürliche Zahlen ver- 
wendet. Eine Beschreibung des Interncodes der Z 22 
und der Eingabe symbolisch adressierter Programme 
mittels eines Adressierprogramms schließen das Buch 
ab. Ein allgemein gehaltener Abschnitt erwähnt die 
Möglichkeiten der Automatisierung der Program- 
mierung (interpretierende und compilierende Technik), 
wobei auf Einzelheiten nicht eingegangen wird. 

Im ganzen erfordert das Buch beim Durcharbeiten 
infolge seiner knappen Schreibweise intensives Mit- 
denken. Es wird den Leser jedoch dadurch anregen, 
die gegebenen Möglichkeiten auszunutzen und selbst 
den Versuch zur Programmierung ihn interessierender 
Probleme zu unternehmen. Die Beispiele stammen 
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alle aus dem Gebiet der klassischen prakti 


matik. Hier wäre eine Auswahl auch aus‘ 
Gebieten der modernen Rechentechnik sehr 
grüßen. Hinweise auf ausführliche Le 


den in einem sorgfältig zusammengestellten Literatur- 


verzeichnis gegeben. 
Berlin K.-H. BAcHMANN 


Handbuch der Physik — Encyclopedia of 
Physies. Band X: Struktur der Flüssigkeiten. 
Hrsg. von S. Flügge. VI + 320 S. m. 41 Abb. Ber- 
lin/Göttingen/Heidelberg 1960. Springer Verlag. 
Preis geb. DM 96,—. 


Dieser Band umfaßt drei Artikel in englischer ne 
che: 1.) Die Struktur von Flüssigkeiten, von Dr. 

BERT $. GREEN (133 Seiten). Nach einer einleitenden 
Diskussion der prinzipiell für Flüssigkeiten zur Ver- 
fügung stehenden Beschreibungsmöglichkeiten werden 
die verschiedenen Modelle diskutiert: Das Gitter- 
modell, das Clustermodell, das Löchermodell und das 
kinetische Modell. Das erste dieser Modelle erinnert 
an die Beschreibung von Festkörpern, das letzte an 
die von Gasen. Die übrigen Modelle berücksichtigen 
die bei Flüssigkeiten besonders auftretenden Schwan- 
kungserscheinungen. Diese Möglichkeiten modell- 
mäßiger Beschreibung werden in den nachfolgenden 
Kapiteln zur Beschreibung der verschiedenen Flüssig- 
keitstypen verwendet, gleichförmige, nichtgleich- 
förmige, quantenhafte Flüssigkeiten. 


2.) Molekulare Theorie der Oberflächenspannung 
in Flüssigkeiten von Dr. Syu Oxo (146 Seiten). In 
diesem Beitrag wird eine Beschreibung der Ober- 
flächenspannung und Adsorption in reinen Flüssig- 
keiten und Lösungen vom molekularen Standpunkt 
aus gegeben. Zunächst werden die phänomenologi- 
schen Theorien der Oberflächenspannung vom Stand- 
punkt der Thermodynamik, der Hydrostatik, der 
Quasi-Thermodynamik besprochen. Ebene und sphä- 
rische (bei Druckdifferenz) Flächen werden gesondert 
untersucht. Dann folgt die eigentliche Darstellung 
vom Standpunkt der statistischen Mechanik. Es fol- 
gen numerische Berechnungen, besonders im Rahmen 
der Gittertheorie. 


Den Abschluß des Buches bildet 3.) die Theorie der 
Kapillarität von Dr. Frank P. Burr (24 Seiten). Die 
bereits von LArLAcE gegebene phänomenologisch 
richtige Beschreibung der Kapillarität wurde von 
GmBBS thermodynamisch fundiert. Die Untersuchun- 
gen vom molekularen Standpunkt aus sind dagegen 
neuesten Datums. 


Dresden W. MackE 


Annual Reviewin Automatic Programming. 
Herausg. R. Goodman, M.A., B. Sc. XI -+ 320 S. 
Oxford/London/New York/Paris 1960. Pergamon 
Press. Preis geb. 63 s.net. 


Der vorliegende Band gibt mit seinen 18 Beiträgen 
Einblick in die Arbeitskonferenz über automatische 
Programmierung digitaler Rechengeräte, die von der 
mathematischen Abteilung der Technischen Hoch- 
schule in Brighton vom 1.—3. April 1959 organisiert 
wurde. Es sollte damit ein erster Rückblick über die 
auf diesem Gebiet speziell in England geleisteten Ar- 
beiten gewonnen und ein Gedankenaustausch über die 
anliegenden Probleme und der dabei aufgetretenen 
Schwierigkeiten ermöglicht werden. Als direktes Er- 
gebnis ist die Konferenzentschließung zu werten, 
zum 1. Januar 1960 in Brighton ein Informations- 
zentrum über Automatische Programmierung zu 
gründen, dessen Aufgaben in der Bereitstellung einer 
Spezialbücherei, in der Herausgabe eines Jahrbuches 
und in der Organisation von Konferenzen liegen, und 


schaffener en hält dieses t 
rende Material zusammen. So gibt der } 
A. D. Booru (Universität London) in Vi 
seiner Eröffnungsansprache und seines S 
mit einer Auswertung der Konferenz wesent 
Situation. Mit einer gewissen | ” 
Gründe untersucht, mit denen weithin Progras 
rungssysteme iert werden. Im Hinblick auf die 
derzeitige Überlegenheit der USA auf dem Feld der 


daß alle diese Fortsetzungsarbeiten auf einem Funda- 


x 


mentaltheorem basieren, das von dem verstorbenen 


englischen Mathematiker A.M. Turıng 1936/37 ge- 
schaffen wurde, Jahre bevor eine moderne ‚Rechen- 
maschine vorhanden war. Daß es verhältnismäßi; 


lange gedauert hat, bis dieser Grundgedanke, daß jede 


Rechenmaschine, wenn sie nur über eine geeignete 
Minimalzahl von Instruktionen verfügt, jede andere 
Rechenmaschine simulieren kann, wie breit auch deren 
Instruktionsvorrat sei, zu der breiten Flut von Pro- 
grammiersystemen geführt hat, findet seine Erklä- 
rung darin, daß die erste Periode, in der Erbauer und 
Benutzer dieselben Personen waren, abgelöst wurde 
durch eine Periode, diein Erweiterung ihres Benutzer- 
kreises auf Techniker und Wirtschaftler eine Vertraut- 
heit mit der Maschine, ja sogar den mathematischen 
Methoden nicht mehr voraussetzen konnte. Ein Nach- 
druck der Pionierarbeit von TURISG, die sich letzthin 
mit dem HıLBErTschen Entscheidungsproblem befaßt, 
wird im Anhang geboten. 


Jena W. KÄMMERER 


U. Bornitz und A. Schubert, Siebenstellige Log- 
arithmen der Fakultäten der Zahlen 0 bis 
12000. 518. Leipzig 1960. Teubner Verlagsgesell- 
schaft. Preis geb. 10,30 DM. 


Viele Formeln der Wahrscheinlichkeitsrechnung 
und der mathematischen Statistik enthalten Fakul- 
täten n!, deren numerische Auswertung bei großem n 
Mühe macht, so daß man oftmals asymptotische Aus- 
drücke benutzen muß. Dadurch gehen allerdings 
Fehler in die Rechnung ein, die meistens nur ziemlich 
roh abgeschätzt werden können. Diese Schwierig- 
keiten werden vermieden, wenn Tafeln für lgn! zur 
Verfügung stehen. Da die bekannten vielstelligen 
Tafeln von €. F. DEGEN und von F.-J. DuvArTE nur 
noch schwer erhältlich sind, ist das Erscheinen des 
vorliegenden handlichen Werkes sehr zu begrüßen. 
Die Berechnung wurde veranlaßt durch Fragen der 
statistischen Qualitätskontrolle und erfolgte auf dem 
digitalen elektronischen Rechenautomaten D1 des 
Institutes für Maschinelle Rechentechnik an der Tech- 
nischen Hochschule Dresden auf 14 Stellen nach dem 
Komma, so daß die Richtigkeit der siebenten Dezimal- 
stelle garantiert ist. 


Stuttgart G. SCHULZ 


I. I. Eterman, Analogue Computers. IX + 2648. 
m. 123 Abb. u. 38 Tafeln. Oxford/London/New York/ 
Paris 1960. Pergamon Press. Preis geb. 50 s. net. 


Das unter der Redaktion von B. H. VENNING aus 
dem Russischen übersetzte Buch von I. I. ETERMAN 
(Moskau 1957) befaßt sich hauptsächlich mit der An- 
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automatischen Programmierung wird mit Stolz betont, 
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wendung elektronischer Analogrechner vom mathe- 
matischen Standpunkt aus gesehen. Nach einem ein- 
führenden Abschnitt über dynamische Systeme wer- 
den im 2. Abschnitt in verhältnismäßig knapper Form 
die Grundrecheneinheiten mechanischer und vor allem 
elektronischer Integrieranlagen behandelt. Im 3. Ab- 
schnitt folgt die kurze Beschreibung einiger sowjeti- 
scher Serienanlagen (u.a. EMU-5, IPT-5, MPT-9, 
MN-7) sowie die Programmierungstechnik mit einer 
Anzahl von Beispielen. Die durch Rechnung ge- 
wonnenen Lösungen dieser Beispiele sind in Abhän- 
gigkeit von der Rechenzeit tin Tafelform zusammen- 
gestellt. Sie stehen somit dem Leser bei der Durch- 
führung von Testaufgaben an Analogrechnern in 
bequemer Form zur Verfügung. Der 4. Abschnitt ist 
spezielleren mathematischen Anwendungen gewidmet. 
So werden Methoden zur Lösung linearer algebraischer 
Gleichungssysteme, Polynomgleichungen und Rand- 
wertaufgaben bei gewöhnlichen Differentialgleichun- 
gen sowie Differenzenverfahren angegeben. Im letzten 


‚Abschnitt werden Genauigkeitsuntersuchungen für 


die Lösung von Differentialgleichungssystemen an- 
gestellt. - 

Wie aus der kurzen Angabe des Inhaltes hervor- 
geht, stehen mathematische Betrachtungen im Vor- 
dergrund. So ist das Werk allen zu empfehlen, die 
sich für die Lösung mathematisch formulierter Pro- 
blemstellungen mit Hilfe elektronischer Analog- 
rechner interessieren. 


Dresden H. ADLER 


H.E. Salzer, D.C. Shoultz, and E.P. Thompson, 
Tables of Osculatory Integration Coeffi- 
cients. 438. San Diego 1960. Convair (Astronau- 
ties) Division of General Dynamics Corporation. 


Ausgangspunkt für die Berechnung der Tafeln_des 
vorliegenden Heftes ist eine n-Punkt-Integrations- 
formel, die man durch Integration der entsprechenden 
speziellen Hrrmiteschen Interpolationsformel (be- 
nutzt werden dabei Funktionswerte f; und erste Ab- 
leitungen f; an n äquidistanten Stützstellen) erhält. 
Vertafelt wurden für n = 2, 3, 4 und 5 die Koeffizien- 
ten von hfi und A2f, für ©,+gh als untere und 
x, + ph als obere Integrationsgrenze als Funktionen 
von p bei fest gewähltem g auf 10 Dezimalstellen, und 
n—ıl 

2 
falls n ungerade, und qg = 0,5, falls n gerade. 
(Die spezielle Wahl von gqin diesen beiden Fällen er- 
folgte, um Symmetrieeigenschaften der Koeffizienten 
voll ausnutzen zu können,) 

Ferner sind in einer Tabelle für alle Koeffizienten 
die exakten Polynomausdrücke (in Abhängigkeit 
von p) angegeben. Eine weitere Tabelle erleichtert 
die Abschätzung des Restgliedes. Schließlich wird an 
einigen Beispielen die hohe Genauigkeit der benutzten 
Integrationsformeln eindrucksvoll vorgeführt. 


zwar für p = q— "0,0 9 +, wobeig=0, 
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I. Babuska, K. Rektorys, F. Vyeichlo, Mathema- 
tische Elastizitätstheorie der ebenen Pro- 
bleme. VIII + 4788. m. zahlr. Abb. Berlin 1960. 
Akademie-Verlag. Preis geb. 96,— DM. 


Das erste und zweite Kapitel (163 S.) benutzen die 
Verfasser zur axiomatischen Einführung der Grund- 
begriffe — direkt für das ebene Problem — und zur 
Formulierung des ersten Problems (Vorgabe der Rand- 
kräfte) und zweiten Problems (Vorgabe der Rand- 
verschiebungen) der ebenen Elastizität mittels holo- 
morpher Funktionen für einfach und mehrfach zu- 
sammenhängende endliche und unendliche Körper 
mit hinreichend glatten Rändern. Die Fälle der un- 
stetigen Belastung und der Belastung durch Einzel- 
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kräfte werden sorgfältig erfaßt. Die Diskussion von 


gemischten Problemen sowie eine Erörterung für 
Eckrandpunkte ist einer späteren Arbeit vorbehalten. 


Die Bestimmung der holomorphen Funktionen des 


Problems wird im dritten Kapitel (164 S.) hauptsäch- 
lich auf die Lösung der Integralgleichung von Lav- 
RICELLA-SCHERMANN (ähnlich der Integralgleichung 
von MUSCHELISCHWILI) zurückgeführt. Aus ihr 
werden einerseits theoretische Folgerungen gezogen 
(Existenzsätze, SAINT-VENANTsches Prinzip, explizite 
Lösung für Kreis und Halbebene), andererseits ist sie 
Ausgangspunkt für Konvergenzbeweise von Nä- 
herungsverfahren; z. B. Ansätze für die holomorphen 
Funktionen in Gestalt von Polynomen und ab- 
brechenden Laurentreihen (Beispiel: Kreisring), ortho- 
normales Funktionensystem bezüglich des Randes 
(Beispiel: Spannung in einem Kranhaken), Nä- 
herungslösungen der Integralgleichung (Beispiel: 
Halbebene mit Kreisloch), Sukzessive Approximation 
nach ScuwArz (Beispiel: Ebene mit zwei Kreis- 
löchern). 

In Kapitel IV (40 S.) werden Probleme mittels kon- 
former Abbildung auf den Einheitskreis übertragen. 
Kapitel V (102 8.) liefert mathematische Hilfsmittel 
aus der Funktionentheorie und der Theorie der Inte- 
gralgleichungen. 

Es ist den Verfassern ein mathematisch exakter 
Aufbau gelungen, der in Verbindung mit zahlreichen 
Vorbemerkungen, Anmerkungen und Beispielen den 
technisch interessierten Mathematiker und den theore- 
tisch interessierten Ingenieur gleichermaßen befrie- 


 digt. 
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A.D. Smirnov, Tables of Airy Functionsand 
Special Confluent HypergeometrieFunctions. 
(Mathematical Tables Series, Volume 7) VII + 260 8. 
Oxford/London/New York/Paris 1960. Pergamon 
Press. Preis geb. £ 5 net. 


Die speziellen Funktionen U;(s, 1), Ux(s,«) und 
V (3,2) (k=1;2), die in dem vorliegenden Band 
(einer Übersetzung aus dem Russischen) samt ihren 
ersten Ableitungen nach s vertafelt sind, sind linear 
unabhängige Lösungen der drei Differentialgleichun- 
gen (DGI.): 

U'(s) +sU(s)=0 


U'(s) + s® U(s) = 0 
Vs) + 5 ( 2) 1 


2 
Ss” 


(1), 
(2) 


2 3 Fa). 0) 


Viele physikalisch-technische Probleme lassen sich 
auf die Lösung einer dieser DGln. zurückführen. 
Ferner können die tabulierten Funktionen dazu 
dienen, die asymptotischen Lösungen der folgenden 
DGl. zu ermitteln: 


|+ era + aly=0 


( 


22 KB dy 
x 


de 


; (A) 


(für ka)>0, 


falls man die ‚Methode der Standard-Gleichungen“ 
von A. A. Doropnıtsyn benutzt. Diese wird an drei 
Typen der DGl. (A) erläutert, deren Lösungen sich 
durch die Lösungen der DGln. (1) bis (3) ausdrücken 
lassen. Für alle tabulierten Funktionen sind Reihen- 
entwicklungen und z. T. Beziehungen zu BESsEL- 
Funktionen mit nicht-ganzzahligem Parameter an- 
gegeben. Einige Diagramme geben einen Überblick 
über den Verlauf der Funktionen. 


Der Tabellenteil enthält in drei Abschnitten: 
1.) Uxls, 1) Ur (sl) 2furzs——6(0.0D10 


(5 Dezimalen bzw. 5 geltende Ziffern, 2. Diff.); 
ferner gewisse Integrale, in deren Integranden die 


r@)>0 mn -aszr=b), 


und 


£3 
« 


Quadrate bzw. das Produkt von U, und D, auftreten: 
Yı1l& 1), LZACH 1), V2(8; 1) für s= 0(0,01) 10 (5 Dez.) 


2.) Urls,a) und U;(s,«) für s = 0(0,01)6 und 
1 1 1 Beleg 
entire 


ARTEN 


3.) Vrls, 9) und V;(s, P,) für s = 0(0,01) 10 
) Vx(8, Po) | \ nn = OLOATE 
(4 Dez.) 


sowie die ersten 5 Nullstellen von V, und Y, für die 
gleichen Werte von »,- 
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M.R.Redwood, (Lecturerin Electrical Engineering. 
Queen Mary College, Univ. of London), Mechanica- 
Waveguides. The propagation of acoustic and ultral 
sonic waves in fluids and solids with boundaries. 
VIII + 300 S. Oxford/London/New York/Paris 1960, 
Pergamon Press. Preis geb. 50 s.net. 


Die Ausbreitung von akustischen Wellen und Ultra- 
schallwellen wird unter den verschiedensten Rand- 
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Ba, } 
bedingungen untersucht. Im ersten Kapi 
eigentliche Theorie aus den Grundpı 
Ten Dr eig 
spezieller. Ran gungen. Zv 
Kapitel über Reflexion und re ee ind 
nachfolgenden Kapiteln fortlaufende Wellen in f} 
gen und festen Platten und Zylindern unters 


Daran schließt sich die Darstellung von Näherungs- 


theorien für kompliziertere Probleme, die Be lur 


von Resonatoren und die Ausbreitung in anisotropen Br. 


Medien. 


Das Buch ist sehr verständlich rn ng Aller- 
dings fällt auf, daß die Vektoranalysis kaum Verwen- 
dung findet und allenfalls nur im früheren Sinne her- 


angezogen wird, um Komponentengleichungen zu- 
sammenzufassen, nicht aber zur Entwicklung physi- 
kalischer Gedankengänge selbst. Das hat leider zur 
Folge, daß die physikalisch wichtigen Ausführungen 
immer wieder durch elementargeometrische Betrach- 
tungen unterbrochen werden. 


Bei den einzelnen Problemen sind Theorie und 


Rechnung vorangestellt, dann folgen Diskussion der 
Ergebnisse und Vergleich mit imenten. Sehr 
wertvoll sind die sehr ausführlichen Literaturhinweise 
hinter den einzelnen Kapiteln. Eine Formelzusam- 
menstellung am Schluß des Buches verbessert den 
allgemeinen Überblick. 
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HELMUT WINKLER 


Elektronische Analogieanlagen 


(Elektronisches Rechnen und Regeln, Band 2) 


1961. X, 242 Seiten — 172 Abbildungen — gr. 8° — Lederin DM 29,— 


In dem Band geht der Verfasser auf das Wesen und die Bedeutung von Analogien und Modellen 
ein und erläutert die Problemerfassung mittels Wirkschaltbilder, Blockschaltbilder, Struktur- 
bilder und Koppelpläne, um letzthin analogieelektronische Lösungswege beschreiten zu kön- 
nen. Er gibt eine Einteilung der elektrischen Analogieanlagen, eine Zusammenstellung aller 
linearen und nichtlinearen Analogieelemente und beschreibt ihr Zusammenwirken. Ferner be- 
handelt er systematisch und ausführlich die mannigfachen Anwendungen elektronischer 
Analogieanlagen mit dem Ziel, den Leser zum Einsatz analogieelektronischer Methoden auf 
seinem Spezialgebiet anzuregen, um den Fortschritt zu beschleunigen, Zeit und kostspielige 
Versuchsaufbauten zu sparen. 

In dem zweiten Teil des Buches schildert der Autor die physikalische und technische Seite, 
die Entwicklungs- und Konstruktionsmerkmale der elektronischen Analogieanlagen, ihre 
Analogieelemente und Baugruppen: Operationsverstärker und elektrische Netzwerke als 
lineare Analogieelemente, nichtlineare Analogieelemente und Funktionsempfänger, Hilfs- und 
Zusatzeinrichtungen wie Meßwandler für den Simulationsbetrieb und schließlich Typen- 
beispiele kompletter elektronischer Analogieanlagen. Abschließend vermittelt er einen histori- 
schen Überblick. Besonders hervorzuheben ist der im Anhang gegebene sehr umfangreiche 


Literaturnachweis, der das einführende Werk wertvoll ergänzt. 
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